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Abstract 

We investigate the asymptotic behavior for type II Hermite-Pade approximation 
to two functions, where each function has two branch points and the pairs of branch 
points are separated. We give a classification of the cases such that the hmiting 
counting measures for the poles of the Hermite-Pade approximants are described 
by an algebraic function h of order 3 and genus 0. This situation gives rise to a 
vector-potential equilibrium problem for measures X, ni, and /U2, and the poles of 
the common denominator are asymptotically distributed like A/2. We also work out 
the strong asymptotics for the corresponding Hermite-Pade approximants by using a 
3x3 Riemann-Hilbert problem that characterizes this Hermite-Pade approximation 
problem. 
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1 Introduction 



1.1 Definition of Hermite-Pade approximants and general state- 
ment of the problem 

Let / = (/i, . . . , /p) be a vector of Laurent series near infinity 

oo „ 

/.■w = E7r' j = h---^p- (1-1) 

k=0 

The Hermite-Pade rational approximants (of type II) 



^ n ^ n 



for the vector / and multi- index n = (rii, . . . , rip) G are defined by 

degP^j <\n\=ni-\ \- Up, 

az)Pn{z) - Q^{z) =: Pg)(z) = (-i^) , z ^ oo, (1.2) 

where the Q^'* are polynomials, for j = 1, . . . ,p. This definition is equivalent to a homo- 
geneous linear system of equations for the coefficients of the polynomial P^. This system 
always has a solution, but the solution is not necessarily unique. In the case of uniqueness 
(up to a multiplicative constant) and in case any non-trivial solution has full degree n, 
the multi-index n is called normal and the polynomial P^i can be normalized as monic 

n 

Pn{z) = Y\iz - Zk,n)- 
k=l 

The Hermite-Pade approximants vr^ provide the best local (near infinity) simultaneous 
rational approximation of the vector (/i, . . . , fp) of Laurent series (11. II) . The construction 
( II. 2p was introduced by Hermite [12] in connection with his proof of the transcendence of 
e. See the papers [531 E21 [13 El IH ESj and the monograph [58j for more details. 

In this paper we study the asymptotic behavior of the diagonal Hermite-Pade ap- 
proximants {n = {n,n)) for two functions /i and /2 with branch points at the points 
Ai = {ai,6i} and A2 = {02,^2} respectively. We say that 

fjeA{C\A,), A, = {a,,b,} (1.3) 

if the Laurent expansion (11. ip is convergent in a neighborhood of infinity and has an 
analytic continuation along any path in C\Aj. A typical example is the function 

fj{z) =\og{{z-aj)/{z-bj)). 

Actually the class of functions we allow is larger and we give a more precise definition in 
Section 2. We shall assume that the pairs of branch points do not coincide, i.e., Ai ^ A2, 
although they might have non-empty intersection. Our goal is 
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• To determine the limiting distribution of the zeros of the common denominator P^, 
which are the poles of the Hermite-Pade approximants: 

^ 2n 

^Pn = l^^^i^ - ^k,n) ^ n ^ OO. (1.4) 



• To obtain asymptotic formulas for the Hermite-Pade polynomials and the func- 
tions P^ . 

• To prove convergence theorems for the Hermite-Pade approximants 

hm = /,(z), zeQ„j = l,2, 

and to describe the domains of convergence Qi and Q2, depending on the location 
of the points oi, 02, &i, &2 G C. 

For the situation under consideration {p = 2 and n = {n,n)) we will use the notation 

= jr = P O^P = O^^^ R^P = R^^^ 7 = 12 

and we assume that P^ is monic. The rigorous definition of the classes of functions 
under consideration and the statement of the results of this paper will be presented in the 
next section. In the following subsections of this introduction we give a brief historical 
review of the analytic aspects of the Pade and Hermite-Pade approximants in order to 
introduce some basic notions and problems. Then we conclude the introduction with a 
general description of the results in this paper. 



1.2 Pade approximants (analytic aspect). Motivation for Her- 
mite-Pade analysis 

The special case of Hermite-Pade approximation (11. 2p for p = 1 (i.e., the best local 
rational approximation to one function near infinity) corresponds to Pade approximation. 
When the coefficients of the Laurent series (11. ip are the moments of a positive measure 
II supported on R 

/(^) = E ^ = / = / di^ix), (1.5) 

then the denominators P„ of the Pade approximants are polynomials orthogonal to the 
powers < n — 1 with respect to fi: 

I Pn{x)x'' dfi{x) = 0, k = 0,l,...,n-l. (1.6) 
Jr 

In the general case when we are dealing with general coefficients in (11.11) the orthogonality 
relations become non-Hermitian or complex. 
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The analytic theory of Pade approximants for the real case fll.Sp is based on remarkable 
classical results. These include Markov's theorem [M] on the locally uniform convergence 
outside the convex hull A of the compact support of /x: 

lim7r„(^)= / zeC\A, 

n-»oo J Z — X 

where K is compact, and the theory of Bernstein and Szego [131 [72] the strong asymp- 
totics of orthogonal polynomials satisfying (11. 6p . 

The analytic theory for the complex case, particularly for functions / G A{£\A), where 
y4 is a finite set of points in C, started to be developed not so long ago. Nuttall [59j has put 
forward the important relation between the maximal domain of analyticity for the analytic 
function / and the domain of convergence of the diagonal Pade approximants. The Pade 
approximants, which are single valued rational functions, approximate a holomorphic 
branch of the analytic function in the domain of their convergence. At the same time 
most of the poles of the rational approximants tend to the boundary of the domain of 
convergence and the support of their limiting distribution models the cuts which make 
the function / single valued. Nuttall has conjectured (and proved for some important 
special cases [59l [63l [60l [38] ) that these cuts have a minimal logarithmic capacity among 
all cuts converting the function to a single valued branch. Thus the domain of convergence 
corresponds to the maximal (in the sense of minimal boundary) domain of holomorphicity 
for the analytic function / G ^(C \ A). The complete proof of Nuttall's conjecture (even 
in a more general setting where the set A has capacity 0) was obtained by Stahl. In a 
series of papers [661 he proved: 

• Existence of a domain VL* such that / is holomorphic in Vt* (/ G H{Q*)) and the 
boundary A = (9f2* has the property 



cap A = min cap dfl. 



• The weak limit of the pole counting measure (11.41) 

z/p„ A A, supp A = A, 



and weak asymptotics for the denominators of the Pade approximants 

lim \og\Pr,{z)\ = -V\z), zen*, 



where is the logarithmic potential 

V\z)= [ log^-rfA(t) 
J \z-t\ 

of an extremal measure A minimizing the energy functional 
= f [ \og-^—d\{t)dX{x) = min 

J J F — t| supp/iCA, /i(A)=l 

among all probability measures on A. 
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The extremal measure possesses the equihbrium properties 



iV = const. 

where denotes the normal derivatives on A (which is a finite union of analytic arcs). 
The conditions (II. 7p characterize the measure A and its support A. These results lead to 
Stahl's main convergence theorem 

lim Tiniz) = f{z), z ^ fl* in capacity, 

n— »oo 

where the convergence in capacity is defined in the same manner as convergence in mea- 
sure. 

The generalization of these notions and results from Fade approximation to Hermite- 
Pade approximation is a very difficult and challenging problem. As we will see, the 
geometry of the domains of convergence and the extremal compact sets where the poles 
of the Hermite-Pade approximants accumulate is much more diverse and complicated. 
The analytic techniques capable of proving the asymptotics and convergence results re- 
quire a significant development in comparison with the methods appropriate for Pade 
approximation and orthogonal polynomials. These circumstances give a good motivation 
and direction for the development of the analytic aspect of Hermite-Pade approximation. 

1.3 Short survey of asymptotic results for general classes of Her- 
mite-Pade approximants 

Perhaps one of the first results on the asymptotics of Hermite-Pade polynomials was ob- 
tained by Kalyagin [H]. For a special class of multiple orthogonal polynomials P^'^''^{x) = 
x^" + • • • generalizing the Jacobi polynomials: 

/ P:;^'^x)x''w{x)dx = 0, fc = 0,l,...,n-l, j = l,2, 

with Ai = [—1, 0], A2 = [0, 1] and weight function 

w{x) = \x + - if, X G Aj, j = 1, 2, 

and for the corresponding functions of the second kind 

R^\z) = ^ ^^^^^w{x)dx, J = 1,2, (1.8) 

he proved the strong asymptotics (Szego type asymptotics) as n — > 00: 

P:i'\z) = Co-"<fo"(^)[Fo(^) + o(l)], ^Gl]o = C\(AiUA2), 
R^\z) = Cr"$T"(^)[F,(^)+o(l)], zen, = C\A„ j = l,2, 
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and the corresponding formulas on the intervals (j = 1,2), where the convergence is 
uniform on compact subsets of the indicated domains. Here 



£ = 0,1,2, 



depend on a, /?, 7 and they are analogs of the Szego function. The main terms $0, '^'i, '^'2 
of the asymptotics are the single valued branches of an algebraic function: 



where the branch points are a,t A = { — 1,0, 1}. We note that $ is independent of 7 
and depends on the supports Ai, A2 of the weight, i.e., on the set A. The function $ is 
a rational function on the three-sheeted Riemann surface 



obtained by glueing the sheets = C \ Aj (j = 1, 2) to the sheet fHo = C \ (Ai U A2) 
so that the upper and lower sides of the cuts on two neighboring sheets are identified. 
The function $ is defined (up to a multiplicative constant) by its divisor (set of poles and 
zeros) , 



and the normalization is chosen so that C0C1C2 = 1. The Riemann surface 9\ and 
rational functions on it play an important role for the asymptotic analysis of Hermite- 
Pade approximants. 

Another important notion of vector potential equilibrium was introduced by Gonchar 
and Rakhmanov in [39j (see also [iQl IH]). Let {Ai, . . . , Ap} be a collection of compact 
sets in C and let D = {di^jY^ be a real symmetric nonsingular positive definite matrix. 
An additional condition on D to be compatible with (Ai, . . . , Ap) is that rfj j > whenever 
Aj n Aj 7^ 0. For a vector of measures 

/!= (/ii,...,/Xp), supp/xjCAj-, j = l,...,p, 

the energy functional I{fl) is defined as 




yK{A) = U 9^1 U 



(1.9) 




(1.10) 



p 



p 




(1.11) 



i=i j=i 



where I{fj,i,fj,j) is the mutual energy of two scalar measures 




The extremal vector measure A, minimizing the energy functional ( II. lip among all fl 
where all fij are probability measures possesses the equilibrium properties 
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Here the vector = {U^, . . . , Up) is called the vector potential of the vector valued 

measure A with respect to the interaction matrix D. 

In the paper [39] Gonchar and Rakhmanov investigated the Hermite-Pade approxi- 
mants f ll.2p for the system of Markov-type functions f ll.Sp 



fj{z) = fij{z) 



X 



fi'j > on Aj C 



where Aj {j = 1, . . . ,p) are non-overlapping intervals 

Ai n Aj= 0, t ^ J, 



:i.l3) 



(1.14) 



where A denotes the interior of the interval A. They proved the weak asymptotics for 
the common denominator where n = {n, . . . ,n) of the Hermite-Pade approximants as 
n — >• cxd: 



1 ' 



-loglPrjW 

n 



z e 



c\J2a; 



3=1 3=1 

where {j = 1, . . . ,p) are the components of the extremal (equilibrium) vector measure 
with matrix of interaction 



d 



'3,3 



di, 



1, 



:i.i5) 



The potentials of the components of the equilibrium measure (after normalization) can 
be harmonically continued through the intervals A* forming the Riemann surface (as in 
(ll.9p where the index j runs from 1 to p and the cuts on the sheets join the endpoints of 
the supports A*). This fact had been noticed in [5]. Thus the notion of rational function 
(ll.lUp on and vector equilibrium problem fll.l2p are equivalent and they are related by 



exp [-V^'{z)) 



C3 



J = l, 




$n(2) 



l! 



where the normalization constants Ci,...,Cp and the equilibrium constants k 
are connected by a linear system of equations. The following convergence theorem was 
proved in [39j: 



lim 




dfij (x) 



z — X 



where 



z e n*, 

ze{C\A))\n*, 

a* = {;.:|$,(^)|>|<|.o(;.)|}. 



j = 1, 



Note that, in view of fll.lOp . fi* 7^ 0. The existence of the non-empty domain C \ fi* 
depends on the input geometry, i.e., on the size and the location of the A^ (j = 1, . . . ,p). 

Thus the results in [SH] show that there are two new phenomena for the asymptotic 
behavior of Hermite-Pade approximants as compared to Pade approximants: 
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1. The components A* of the support of the pole counting measure do not correspond 
to cuts making the functions fj{j = l,...,p) holomorphic. We call this the pushing 
effect: it might happen that \ A* 7^ for some j e {1, . . . 

2. The appearance of domains of divergence inside the domain of holomorphicity of fj-. 
it might happen that (C \ Aj) \ fi* 7^ for some j G {1, . . . 

These two phenomena are related by 

A,\A*^0 J (C\A,.)\^,V0- 

The system f ll.l3l) - fll.l4l) was introduced in 1919 by Angelesco [2j as a system for which 
all the multi-indices of the Hermite-Pade approximants are normal, and this system was 
later rediscovered in [55J. Another system of Markov- type functions (11. 5p with normal 
diagonal multi-indices for the Hermite-Pade approximants was introduced by Nikishin in 
[56]. A system (11.131) is a Nikishin system of order p if Aj = A for j = 1, . . . ,p and 
dnj/dfii [j = 2, . . . ,p) have analytic continuation from A and form a Nikishin system of 
order p — 1 with respect to another interval F for which F n A = 0. The asymptotic 
behavior of the denominators of the Hermite-Pade approximants for a Nikishin system is 
similar to the behavior of Pade approximants in the sense that the two phenomena for 
Angelesco systems do not appear (see [56l[571[20l[3ll[68]). However, for a Nikishin system 
a new effect appears for the functions of the second kind (11.81) 

E!i\z) = Pn{z)f,{z) - (^{Z), J = 2, . . . ,p. (1.16) 

They have extra zeros which accumulate on the interval F and are dense on this interval 
as 00. These are extra interpolation points for the Hermite-Pade approximants, 

apart from the interpolation condition at 00. 

To conclude this survey of the results for the real case, i.e., Hermite-Pade approxi- 
mants for a vector of Markov type functions, we mention the recent paper [5T] on mixed 
Angelesco-Nikishin systems defined by a graph-tree, and the papers [J, 5J where the strong 
Szego-type asymptotics of the Hermite-Pade polynomials for Angelesco systems and Nik- 
ishin systems was obtained. 

The analytic theory of Hermite-Pade approximants for the complex case has been 
initiated by Nuttall. In the two pioneering papers [121 EI] of 1981 he obtained some 
asymptotic formulas for Hermite-Pade approximants to functions with separated complex 
branch points [12] (a complex analog of an Angelesco system) and to functions meromor- 
phic on the same Riemann surface [61] (i.e., functions with the same set of branch points, 
like a Nikishin system for the real case). The results of [12] were verified by some heuristic 
considerations and numerical experiments, and the paper [HI] contains rigorous theorems. 
In his fundamental work [62j of 1984, Nuttall made an attempt to formulate a general 
conjecture about the asymptotic behavior, as n — > co, of the diagonal n = {n,n, . . . ,n) 
Hermite-Pade polynomials. On the basis of his conjecture lies a (p + l)-sheeted Riemann 
surface like (11.91) ) which depends on the set of p functions which are being approximated. 
He showed how to determine this Riemann surface for some special classes of functions, 
but the general case was left as an open problem. Nevertheless, assuming the existence 
of the appropriate Riemann surface fH, he conjectured that the strong asymptotics can 
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be described by solutions of some boundary value problem on For the main term of 
the asymptotics one would have 



\Pn{z)\^^'" ^\'^~\z)\=exp{~Re G{z)) , n oo, (1.17) 

where G is an Abelian integral of the third kind with logarithmic poles at oo'-^-', (£ = 
0, 1, . . . ,p) with residues 

f-plog^ + 0(l). .^oo(»). 

\log2 + 0(l), Z^C0''K J = l,...,p, 

and elsewhere G is analytic in the local variable. If the genus of is greater than zero, 
then an additional condition on G is imposed: all the periods of G are purely imaginary. 
Such a function is unique up to an additive constant and Re G{z) is a single valued 
function on 91 (see, for example, [68]). We note that the condition genus(fH) = implies 
single- valuedness of $ in f ll.l7p . therefore $ is a rational function on IK, uniquely (up to 
a multiplicative normalization) defined on by its divisor f ll.lOp . 

After Nuttall's results there were practically no other rigorous results for the complex 
case of Hermite-Pade approximant^. One of the reasons was the absence of suitable 
techniques for the analysis of strong asymptotics of non-Hermitian (complex) orthogonal 
polynomials which could be adopted to Hermite-Pade approximation. There was recently 
substantial progress in proving new results for the strong asymptotics of orthogonal poly- 
nomials by means of a matrix-valued Riemann-Hilbert method. The method is based on 
the reformulation of the definition (II. 6p of the orthogonal polynomials in terms of a 2 x 2 
matrix- valued Riemann-Hilbert problem (due to Fokas, Its, and Kitaev [371 [36]) ^^id the 
steepest descent analysis of this Riemann-Hilbert problem for n — oo (due to Deift and 
Zhou [31]). This method was initially designed to study asymptotics for integrable PDEs 
and was later applied to prove asymptotic results for polynomials orthogonal on the real 
axis with respect to real valued analytic weights, including varying weights (depending 
on n) [m [26|, |271 [281 [291 [30l [IHl [U] and related questions from random matrix theory. It 
has later been noticed [6l [Til 05] that the method also works for the non-Hermitian or- 
thogonality in the complex plane with respect to complex valued weights. In [H] multiple 
orthogonality for Hermite-Pade polynomials was reformulated in terms of a (p+1) x [p+ 1) 
matrix-valued Riemann-Hilbert problem. 

In this survey subsection we introduced the main players for the asymptotics of the 
complex case of Hermite-Pade approximation in the historical order of their appearance 
on the scene. They are 

• An appropriate {p + l)-sheeted Riemann surface 9^. 

• Vector potential equilibrium problems. 

• A standard Abelian integral fll.181) and rational functions on 9^. 

• Analysis of a (p+ 1) x (p + 1) matrix-valued Riemann-Hilbert problem as a method 
for proving the asymptotic results. 



^except for the special case of Hermite-Pade approximation to the [HI [33l [35l HH |49l [50l [Ml EOl 

iniiiais]. 
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Concluding this survey of results which have important and direct influence on this 
paper, we would like to mention the papers [6], [3 HUl [16], |l7l [251 [SOI [SH [52] which have 
been written during the work on the present paper and in which some of the ideas and 
methods elaborated here have already been implemented. 



1.4 General description of the limiting behavior of Hermite- 
Pade approximants 

The analysis of numerical computations of the zeros of the polynomials P^, which are 
the common denominators for the Hermite-Pade approximants to functions of the class 
( II. 3p . shows that the support for the limiting zero distribution has a different geometry, 
depending on the position of the branch points. In Figures I1.1H1.5I we present the results 
of the computations for the functions 



fj{z) = log 



J = 1,2. 



There are several typical patterns for the support. We describe those that are related to 
a Riemann surface of genus zero. 

Case I. When the pairs {ai, bi} and {02, ^2} are 'far away' from each other, the zeros of 
Pfi accumulate on two disjoint arcs Ai and A2, which are joining the branch points 
Oi and 61 for the function /i and the branch points 02 and 62 of the function /2, 
respectively (Figure [LTl) . Each of the arcs Ai and A2 accumulates half of the zeros 

of Pn. 




Figure 1.1: ai = 0, 61 = i, 02 = 2, 62 = (3 + 3i)/2 (case I). The zeros of P„^„ accumulate 
on two disjoint arcs Ai and A2 that connect the branch points. The zeros of P2o,2o are 
indicated by o. 
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Case II. When the pairs {ci, bi} and {02, 62} are closer to each other, then it may happen 
that the zeros of accumulate on a set A that connects all four branch points ai, 
bi, 02 and 62 as shown in Figure [L2l In this case the support can not be split into 
two separate pieces of equal mass. 




-0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 



Figure 1.2: ai = z, 61 = 0, 02 = 1/2 + z, 62 = 1/2, (case II). The zeros of Pn,n accumulate 
on a contour A that connects all four branch points. The zeros of ^20,20 are indicated by 



Case III. For certain configurations of pairs {ai,bi} and {02? ^2}; the zeros of accu- 
mulate on a set Aq that consists of two disjoint arcs A J and A2 (as in case I), but 
contrary to the case I, the arcs are not joining all branch points, see Figure [L3l One 
of the branch points (61 in Figure fTISi) does not belong to Aq. The arc Al connects 
Oi with a point b* which is different from bi. The other arc A2 connects 02 with 62- 
Both arcs accumulate half of the zeros of P^. 



0.5 



0.0 



-0.5 



1 1 1 1 1 






02 




b* 


r 








»u u n 1 n ^ e e a a- 






ai 


b\ \ 








b2 


1 1 1 1 1 



-0.5 0.0 0.5 1.0 1.5 2.0 2.5 



Figure 1.3: ai = 0, 61 = 15/8, 02 = 2 + i/4, 62 = 2 — i/8 (case III). The zeros of P„_„ 
accumulate on two disjoint arcs A^ and A2. The branch point bi is not contained in A^. 
The zeros of P2o,2o are indicated by o. 
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The above three cases concern situations where the four branch points aj, bj are all 
distinct. The final two cases deal with situations where {ai,6i} and {025^2} have a 
common point, and we take it so that 



bi = b2 = b. 

Case IV. The pairs {ai, b} and {02, b} are such that the zeros of Pa accumulate on a set 
A that connnects all three branch points. 




-1.2 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 
Figure 1.4: ai = —1, 6 = 0, 02 = i/2 (case IV) 



Case V. Certain positions of the pairs {ai, b} and {a2, b} are such that the zeros of Pa 
accumulate on a set Aq that consists of two disjoint arcs A* and A2 as in case III. 
Thus, A2 connects 02 with b, but the other arc A* connects ai with a point b* which 
is different from b, and which does not belong to A2. Some of this situations can 
be realized as limiting cases of case III where the point bi tends to 62- These cases 
are not part of the case V, but rather belong to case III. The case V contains the 
situations that cannot be realized as limiting cases of case III. 




-1.2 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 



Figure 1.5: ai = —1, 6 = 0, 02 = 0.32i (case V) 

There are also 'higher genus cases' such as the one shown in Figure [L6l Here the zeros 
of Pfi accumulate on a set A consisting of three disjoint arcs. Each of the branch points 
is contained in A, but only two of them (in the figure it is 02 and 62) are on the same arc 
A2. This arc accumulates half of the zeros of P^. The other half are on the two remaining 
arcs. We will not treat the higher genus cases in this paper. 
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1.2 



ai 



0.6 



02 



0.0- 



-0.6 



-1.2 



bi 



T 



-0.2 0.0 
Figure 1.6: ai = i, bi 



0.2 0.4 0.6 0.8 1.0 1.2 

— 1, a2 = 1 + i/2, 62 = 1 — i/2 (higher genus case). 



There are also critical cases when there is a transition from one case to another case, 
e.g., the critical case where the two arcs Ai and A2 have one common point, which is the 
transition from case I to case II. 

The cases I and III present a complex generalization of an Angelesco system f ll.lSp - 
(I1.14p . As we will show, the limiting distribution of the poles of the Hermite-Pade approx- 
imants for these cases is the equilibrium measure of the vector potential extremal problem 
fll.lip - fll.12l) with the Angelesco interaction matrix (11.151) . and the support A := Ai U A2 
of this measure is characterized by a vector analog of Stahl's symmetry property 



dul 



dul 



on Aj, j 



1,2. 



The cases II, IV and V have features of a Nikishin system because of the presence of the 
common piece and we cannot separate A into two disjoint pieces. For these cases a new 
curve E appears on which the finite zeros of the function of the second kind (i.e., extra 
interpolation points) accumulate and an extremal vector potential problem describes both 
the limiting distribution of the poles of the Hermite-Pade approximants and the extra 
interpolation points. 

The three-sheeted Riemann surface which is appropriate for the cases I, III and V 
has four branch points at the end points of the components of the support of the vector 
equilibrium measure (for case I they are the branch points of the functions which are being 
approximated). Thus this Riemann surface has genus zero. The three-sheeted Riemann 
surfaces for the cases II and IV have a more complicated sheet structure but their genus 
is also zero. 
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1.5 Objectives, structure, main results, and tools of the paper 



In this paper we prove the asymptotic formulas and the convergence theorem of the Her- 
mite-Pade approximants for the cases when the appropriate Riemann surface is of genus 
zero. In a forthcoming paper we plan to describe the cases when the appropriate Riemann 
surfaces are of genus one and two. In this paper we will handle cases I, II, III, IV, V that 
give rise to a genus zero Riemann surface as described above. 

Rigorous definitions of the class of functions which we are approximating corresponding 
to the geometrical cases and statements of the results will be presented in the next section. 
The rest of the paper is devoted to the proofs of the theorems. Here we present a general 
description of the obtained results. 



1.5.1 Geometry of the problem 

Our approach to the geometry of the problem, i.e., finding the appropriate Riemann 
surface, classification of the geometrical cases, etc., is based on the algebraic function h 
defined by 

where the polynomial 114 is defined by the input parameters (the branch points of the 
functions we are approximating) 

n4(z) = {z- ai){z - bi){z - a2){z - 62), 

and the polynomials P2 and Pi contain two parameters 

P2iz) = z^ - diz + d2, 

P^[z) = z-c, (1.20) 

2c + fli + 61 + 02 + &2 

"1 = 1 

3 

which can be determined by the input data ai,a2,bi,b2, using the information about 
the geometrical case to which the input data belongs. The function h has three branches 
ho, hi, h2 which we fix at infinity as 

2 

hoiz) = — + ••• 

z^oo. (1.21) 

hj{z) = - + ■■■, J = 1,2. 

The function /i is a rational function on its Riemann surface 9^. It has 

• four poles at the points on for which the projection vr : 9i ^ C are the points 
{01,02,61,62}: 

/i(0 = oo ee^R:7r(0 G {01,02,61,62}; (1-22) 

• three zeros at infinity on each sheet and one zero at the point for which the projection 
is the parameter c from (11.201) : 

^ 1^--'' (1.23) 
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If the zero at c cancels a pole from {oi, 02, 61, &2}; for example when c = 61, then the 
equation (11.191) for the algebraic function h reduces to 

^ {z - ai){z - a2){z - ~^ {z - ai) {z - {z - ^' ^ ^ 

where 

, ai + 02 + &2 

^= 3 • 

The discriminant V of the equation fll.24p is 

^ _ {z - d)^ - {z - ai){z - a2){z - 62) 
{z - ai){z - a2){z - b2) 

hence the function h from fll.24p has branch points at ai, 02? &2 (the poles of h) and at the 
point 

b* = ^^^2^^ - (1.25) 

aia2 + 0162 + 0.2^2 - 3(i2 

The Riemann surface corresponding to equation (11.241) has genus zero. This Riemann 

surface was considered for the first time by Nuttall [621 [12] by means of another equation, 

and independently by Kalyagin [8]. In our work this Riemann surface will be an 

appropriate 9^ for the geometrical cases III and V (for Case V the functions we are 

approximating have a joint branch point, i.e., bi = b2 = b in (II. 3p ). 

For the geometrical Case I the appropriate Riemann surface is defined by the general 

equation (I1.19P where the two unknown parameters in (11.201) are determined by the two 

conditions 

1. the genus of IK is zero, 

2. the monodromy condition 

ho e H(C \ (Ai U A2)), hj e H(C \ Aj), J = 1, 2, Ai n A2 = 0. 

Here Aj is an arbitrary Jordan arc joining the points aj and bj {j = 1,2). The discriminant 
V of the equation (I1.19P is 

114 

where the polynomial T> has degree 4 due to (11.201) . Condition 1 above implies that 

V{z) = const {z - zif{z - Z2f, (1.26) 

which gives a system of two algebraic equations (of high order) for the determination of 
the two unknown parameters in (ll.20p . This system of algebraic equations has several 
solutions with different monodromy properties, for example, G H{C \ Ai), hi E H{C \ 
(Ai U A2)), h2 G H{C \ Ai), hence condition 2 above chooses the right solution. We will 
not work with this system of algebraic equations because it is too cumbersome. Instead we 
introduce in the next section a substitution for (ll.20p which automatically fulfills (11.261) . 

Thus the geometrical part of the problem consists of the determination of the algebraic 
function h starting from the input branch points Aj = {aj, bj}, [j = 1,2). As we already 
mentioned above, we characterize in this paper the position of Ai and A2 which guarantee 
that the algebraic curve (I1.19P is of genus zero. 
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1.5.2 Standard functions for the asymptotics 

Since ttie geometrical analysis gives us the algebraic function h and its Riemann surface 
we can define the functions which allow us to state the asymptotical results. The 
Abelian integral, see (11.181) . is defined as 

G{0= [^h{z)dz, (1.27) 

and the function 

$(0=expG(0, (1.28) 

is a single valued (rational) function on 9^, which is a consequence of the fact that the 
genus of £H is zero. The local selection of the branches of the algebraic function h in (11.211) 
gives us the local definition of the branches for the algebraic function $ (see (ll.lOp ): 

1 z 
$0(2;) = + ■ ■ ■ , <l>j{z) = — + ■■■ , 3 = 1,2, z^oo, 

C1C2C0 = 1, Ci > 0. 
Using the union of the analytic curves defined by 

r = {zeC: \<^e{z)\ = for some < £ < k <2}, 

we define the holomorphic branches of $ (and respectively h) globally: 

Gi/(C\7£), 7iCT, i = 0,1,2. 

Details will be given in the next section. 

Thus the appropriate Riemann surface = IHo U U 9^2 can be realized as three 
sheets of the extended complex plane cut along the contours •jf 

and pasted through 

WenWk = je,k, ^^k, = 0,1,2, 
so that the upper and the lower sides of the cuts on two neighboring sheets are identified. 

1.5.3 Asymptotic results and convergence 

Using the global definition of the branches of the standard algebraic functions h and $ 
(depending on the geometrical case) we can sketch our asymptotic results. The limiting 
distribution (ll.4p of the zeros of the diagonal n = {n, n) Hermite-Pade denominator 
and a strong asymptotic formula can be written in a unique way for all geometrical cases 
under consideration. In the present paper we prove that 

* 1 

where the real valued measure A, which is of total mass two, is given by 

dKz) = 7^ [iho)+ - (ho)-] dz, z e 70, (1-29) 
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where the subscript + or — as usual denotes the hmiting value of the function taken from 
the left or right, respectively, when traversing the contour according to the orientation 
determined by the complex line element dz. The strong asymptotics for is 

Pr-t{z) = Co-"<f "(^) {F,{z) + 0(l/n)) , 2 G C \ 70, 

where Fq G if(C \ 70) is an analog of the Szego function defined by mean of a certain 
boundary value problem (see details in the next section). 

The asymptotic formulas for the functions of the second kind f ll.Sp and fll.l6p essen- 
tially depend on the geometrical case. For the cases I and III we prove that 

4'^W=C-"$7"(F,(^) + 0(l/n)), ;.gC\7„ J = 1,2, (1.30) 

for certain functions Fj G if(C \ 7^) (j = 1,2). For the case II, IV and V the answer 
for the functions of the second kind is more involved. A general feature for these cases is 
that the asymptotic formula fll.301) remains valid in a neighborhood of 00 which is smaller 
than C \ 7j, j = 1, 2, and around the 7^ (j = 1, 2) there may appear a domain where one 
of the main terms of the asymptotics of {$1, $2} changes from one to the other, so as a 
result the functions of the second kind on the boundary of this domain have an oscillatory 
asymptotic behavior which leads to an accumulation of its zeros there. 
From (11.161) we have that 

^'^^ Pn{z) Priiz) ' 

hence our asymptotic results give a complete picture of the convergence of the Hermite- 
Pade approximants with a description of the possible regions of divergence and the sets 
of the accumulation points of the extra interpolation points. 

1.5.4 Tools 

To prove the asymptotic results we start with a 3 x 3 matrix-valued Riemann-Hilbert 
boundary value problem characterizing Hermite-Pade approximants for two functions, 
which was proposed in fT^. Then, using the information about the geometry of the 
problem, we develop the steepest descent method of Deift and Zhou [3T] which was already 
successfully used for 2x2 matrix-valued Riemann-Hilbert problems. 

We can say that for the geometrical cases I and III the jump matrices in the initial 
Riemann-Hilbert problem have a block structure and therefore most of the steps of the 
asymptotic analysis of the solution of the 3x3 matrix-valued Riemann-Hilbert problem 
can be reduced to the 2x2 problem, which has already been developed. As a new feature 
we like to mention the procedure of finding the explicit solutions for the matrix Riemann- 
Hilbert problem with non- varying jumps (i.e., jumps which do not depend on n), reducing 
the matrix problem to a boundary value problem on the corresponding Riemann surface. 

For the geometrical cases II, IV and V the jump matrices of the Riemann-Hilbert 
problem do not possess this block structure. We introduce a new decomposition of the 
jump matrix to block structure jump matrices (which can be treated by the traditional 
local decomposition) and a jump matrix with exponentially growing non-diagonal terms. 
Nevertheless, due to the analyticity of the solution of the Riemann-Hilbert problem, the 
contour on which this growing jump occurs can be moved in the domain where the growing 
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terms of the jump become exponentially decaying. We call this new decomposition a global 
opening of the lenses. This new procedure brings about new curves which initially have not 
been present in the statement of the problemo For example, we will discover the analytic 
curves on which the extra interpolation points accumulate. We also like to mention a 
wonderful picture: when the contours of the jumps make their global movement they may 
meet each other or one may pass through another and the corresponding jumps interact: 
we can see interference, transparent penetration or even annihilation. 



2 Rigorous definitions and statements of the results 

2.1 Class of functions and reformulation of the Hermite-Pade 
approximation problem 

Now we will be more precise in the definition of the class of functions (11. 3p we are ap- 
proximating and we state a matrix Riemann-Hilbert problem as a reformulation of the 
Hermite-Pade approximation problem. 

Let a and b be points in the complex plane and let A be a Jordan rectifiable arc joining 
a and b: 

a,beC, A = {A(t) : t G [0,1]}, A(0) = a, A(l) = 6. (2.1) 
We will consider functions / of the form 

where w is some "nice" function on A as specialized in the following definition. Note that 

w{0 = U{0-f-{0, eeA, (2.3) 

where f± denote the boundary values of / from the left and right using the orientation 
on A from a and b. 

Definition 2.1. Let a, (3 > —1 and let Qhe a. domain containing the arc A joining a and 
b. Then we say that the function / given in (12. 2p belongs to the class A{a, a;b, P;fl), 

f eA{a,a-b,p-Q), (2.4) 

if the function w given in (12.31) satisfies 

^(0 = ^o(0(«-er(e-&)^ a,P>-l, 
Wo,l/woeH{Q). 

Here we fix a holomorphic branch on A C ^2 of (a — z)°'{z — b)^. Thus a function 
/ G A{{a, a]b, (3; fl) has analytic continuation across the arc A into the domain Q situated 
on the next sheet of the Riemann surface of /. Moreover in the representation (12. 2p of 



^The phenomenon of global opening of lenses was first discovered during the research leading to this 
paper. It turns out to be a general feature for higher order RH problems and was also used in the papers 



mm- 
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f{z) we can choose any arc A in i7 that is homotopic in Q to the original arc, as long as 
we keep the endpoints fixed. 

For functions /i and /2 satisfying (I2.2p - fl2.3l) . the definition of the Hermite-Pade ap- 
proximants ( 11. 2p is equivalent (by Cauchy's theorem) to the determination of the multiple 
orthogonal polynomial Pni,n2 which satisfies 



/ i^„„n2(OeS(Oc?e = 0, A; = 0,1 



' • • • ' '"J 



- 1, j = 1,2, 



(2.6) 



and the functions of the second kind are given by 



^ni,n2(0Wi(0 



27ri Ja^, z 
Therefore the error of the approximation is 

Qni,n2(2^) -^"i,n'2('^) 



J = l,2. 



(2.7) 



-f 711,712 (^) 



711 ,712 (^) 



(2.8) 



We note that the right hand side of (12.81) is independent of the normalization of Pni,n2- 
We will fix the normalization of Pni,n2 by choosing it to be a monic polynomial of degree 
rii + n2. Then 

P„„„,(z) = z"^+"^ + ---, (2.9) 



and 



^ni,7i2 * 
1 



?^7ii,n2 

However, the existence of Pni,„2 satisfying (12.61) and (12. 9p is not guaranteed: we recall 
from the introduction that this is only guaranteed for normal multi-indices. 

For normal multi-indices the Hermite-Pade polynomials -Pni,n2 ^ind the functions of 
the second kind (12. 7p can also be defined by means of a matrix-valued Riemann-Hilbert 
problem [71]. We state the problem for the diagonal multi-indices (the case of our interest) 



(i) 



m 
-1 
2tH 



^ — > OO 



711,712 



(2.10) 



n = {n, n), ni = {n — 1, n), n2 = {n, n — 1). 
If n, ni,n2 are normal indices, then the 3x3 matrix- valued analytic function 



(2.11) 



/ P,{z) Rf{z) Rf{z) \ 



Yiz) 



miPn^{z) miRf^{z) miRl^'{z 
\m,P^,{z) m,R^ll{z) m^Rf^iz) j 



zGC\Ao, Ao = AiUA2, (2.12) 



with mi = m!'}^ and m2 = m)^_^ (see (I2.10p ). is a solution of the following Riemann-Hilbert 



,(2) 
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problem: 



1. F G H^^^{C \ Ao) with boundary values Y± on Ao- 

2. Jump relation on Aq 

3. Normalization near infinity: 



Yiz) 



Here the jump matrix is 



/ + ( - 

z 












oo. 



(2.13) 



1 







1 



(2.14) 



where we assume that 



Wi 

W2 








on A2 \ Ai 
on Ai \ A2, 



and / is the 3x3 identity matrix. If we suppose in addition that the entries of the jump 
matrix (I2.14p belong to the class (12.51) . then the matrix Riemann-Hilbert boundary value 
problem (12.131) . supplemented with conditions at the endpoints of Aq = Ai U A2, has a 
unique solution. The endpoint conditions when 2; — > ai are 



Viz) 



0(1 

V0(i 

/0(i 
0(1 

V0(i 

/0(1 
0(1 
V0(i 



0(|^-ai|°i) 0(1)\ 

0(|2-ai|"O 0(1) : 
0(|2-ai|"i) 0(1)/ 

0(log|^-ai|) 0(1)\ 
0(log|^-ai|) 0(1) 
0(log|z-ai|) 0(1)/ 

0(1) 0(1)\ 

0(1) 0(1) , 
0(1) 0(1)/ 



-1 < ai < 0, 



«! = 0, 



(2.15) 



«! > 0, 



Here oi is assumed to be an endpoint of Ai and not an endpoint of A2. If cii coincides 
with an endpoint of A2 (e.g., ai = 02), then we will assume that 



Oi = 02 



«! = 02, 



(2.16) 
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and the endpoint conditions as ^ — >■ Oi become 



Yiz) 





/0(1) 


0i\z 






0(|z- 


ai 








0(1) 


0{\z 








ai 


















ai 








/0(1) 


0(log 


[\z - ail) 


0(log 




— ai 




< 


0(1) 


0(log 


- ail) 


0(log 




-ail) 




^0(1) 


0(log 


I \z - ail) 


0(log 




— ai 


IV 




/0(1) 


0(1) 


0(1)\ 














0(1) 


0(1) 


0(1) 
















0(1) 


0(iV 













-1 < ai = 0:2 < 0, 



Oil = a2 = 0, 



«i = «2 > 0, 



(2.17) 



When z approaches an endpoint other than oi then the endpoint conditions are similar 
to (l^J[^ - (ETri) . We then have 



Proposition 2.2. Suppose the Riemann-Hilbert boundary value problem (12.131) has the 

matrix jump W given by fl2.14p with entries satisfying (12.51) and endpoint conditions 
fl2.15p or fl2.16l) - fl2.17p (and similar conditions at the other endpoints). If this Riemann- 
Hilbert problem has a solution, then this solution is unique and its entries fl2.12p are the 
Hermite-Pade polynomials (multiple orthogonal polynomials) (12.61) and the functions of 
the second kind (12.71) corresponding to the normal index n = {n, n) of the Hermite-Pade 
approximants for the functions 



fj G A.{aj, 01 j] bj, Pj] Qj^ 



1,2. 



(2.18) 



Proof. See [74] for the relation between this 3x3 matrix- valued Riemann-Hilbert problem 
and multiple orthogonal polynomials. The uniqueness when the endpoints conditions are 
imposed follows in a similar way as in ^7|. □ 



Our goal is to find a solution of the problem (I2.13p - (l2.17p for large n. This would 
give the normality of the Hermite-Pade approximants for large n and their asymptotics 
as n — i> 00. For this purpose we apply a steepest descent method for asymptotic analysis 
(as n 00) of the matrix Riemann-Hilbert problem. One of the main ingredients of this 
method consists in finding a location for the arcs Ai and A2 within the domains Qi and 
Q2 such that the jump matrices in the Riemann-Hilbert problem admit a factorization as 
a product of matrices which tend exponentially fast (as n —>■ 00) to the identity matrix 
outside A = Ai U A2 and a matrix independent of n. The location of these arcs depends 
on the position of the points {ai, a2, &i, &2}- Thus the form of the asymptotics of the 
solution of the Riemann-Hilbert boundary value problem and the method of obtaining it 
depend strongly on the position of the points {ai, 61, a2, 62}- 

We will assume throughout that the domains fli and Q2 are such that they contain 
the "optimal" arcs Ai and A2. We also assume that the domains are large enough so 
that the deformations we need to do in the course of the steepest descent analysis can be 
performed within Qi and ^2- This is in particular important for the cases II, IV, and V 
where one of the steps involves global deformation of contours. 



22 



2.2 Geometry of the problem. Cases I and II 

In this section we define a class of positions of the points 

A:= {ai,bi;a2,b2} (2.19) 
which we call the geometrical cases I and II 

A e lull. 

This class characterizes analytic functions (12.41) 

/, E A{a„a,;b„P„n^), j = 1,2, (2.20) 

whose Hermite-Pade asymptotics are described by means of algebraic functions of the 
third order with the only branch points Oi, hi, 02, &2- We define the two classes I and 
II and we fix for the input data (I2.19P the global branches of the algebraic functions 
needed for the presentation of the asymptotic results and we introduce measures for the 
description of the limiting behavior of the zeros of the Hermite-Pade polynomials. 

2.2.1 New coordinates for the input data 

The algebraic function h (see (11.191) ) associated with the input data (12.190 satisfies the 
equation 

hHz)~3^hiz) + 2^ = 0, (2.21) 
114(2;) 114(2;) 

and plays a key role for the classification of the geometrical cases. We recall from the 
introduction that 114(2;) = (2; — ai ) (2; — 61) (2; — 02) (2; — 62) and the two unknown parameters 
of the monic polynomials Pi and P2, of degree one and two respectively (see (11.201) ). have 
to satisfy — for the cases under consideration in this paper — the condition that the 
genus of h is zero. It follows from this that the discriminant T> of the equation ( I2.2ip 

has zeros of even multiplicity, i.e., (I2.19P are the only branch points of h. As we already 
mentioned in the introduction, this gives us two algebraic equations (of rather high order) 
for the determination of the unknown parameters of Pi and P2. In practice it is better to 
consider the algebraic function h with a special class of coefficients in (12.211) satisfying the 
condition that the zeros of the discriminant of h have multiplicity two. In other words, 
starting from known coefficients in (12.211) . we will obtain the set A in (I2.19P as the set of 
zeros for the polynomial 

I1a{z) = {z- ai){z - bi){z - a2){z - 62). (2.22) 

To realize this, we introduce new coordinates {k,p,s,c) in for the input data. The 
coefficients in (I2.2ip are given by the following 
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Proposition 2.3. Suppose we are given 



5 = (A;,p,s,c) G C^ (2.23) 

and we define 

Pi{z)=z-c, P2iz) = P^{z) + 2pPi{z) + s\ (2.24) 
Then _ 

U,{z) = with Viz) = (kP^iz) + 3psPiiz) + s^f, 

is a polynomial of degree 4 in the variable z: 
n4 = Pi^ + 6pPi^ + (3s2 + 12p2-A;2)p2^(8p3 + i2s2j9-6psfc)Pi + 3/ + 3p^s^-2s^A;. (2.25) 
The proof of Proposition 12.31 is given in Subsection I3.1.1[ 

So given /c,p, s, c we construct the polynomials Pi, P2, and 114 as in fl2.24p and (12.251) 
and we use these polynomials in the definition (12.211) of the algebraic function h. Then 
it follows from Proposition 12.31 that D is a perfect square, and (I2.2ip has branch points 
at the zeros of 114 only, which we denote by ai, 61, a2, 62 (taken in some order) and we 
assume that they are all distinct. This defines a mapping 

B = {k,p,s,c) ^A = {a^,b^,a2,b2} (2.26) 

Thus, if we use the coordinates B in (I2.23P then we can determine both the set A of four 
points in (I2.19P and all the coefficients of the equation (I2.2ip of the algebraic function 
h, which has its branch points in A. Generically, this function h has genus zero, that is, 
with the exception of a lower dimensional subset of the four dimensional {k,p, s, c)-space. 
However, this class also contains functions of higher genus (for example, functions with 
branch points of the third order). These degenerate examples are not connected with our 
task and we exclude them from our considerations. 



2.2.2 The function $ and the contour P 

Since the algebraic function h has genus zero, the exponential of the Abelian integral 



$ = exp / h{z) dz 



(2.27) 



is also an algebraic function with the same Riemann surface as h. The equation fl2.27p 
defines $ up to a multiplicative constant. 

Proposition 2.4. For a suitable normalization we have that the function $ satisfies the 
equation 

^\z) + q^{z)^\z) + q2{z)^{z) + go = 0, 
where qj are polynomials of degree < j {j = 0, 1, 2) given by 



-eVSpPiiz) + ^\{k- 2V3p)K+ + {k + 2V3p) 

9s2 _ 4^2 ^ 3gpS 



-K+K^ ( Pl{z) + ApPi{z) + 

2v^ 
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(2.28) 



■(3s - 2kfKW 



243 
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where Pi{z) = z — c and 

K± = 3s + k± 3VSp. (2.29) 
The proof of Proposition is in Subsection I3.1.1[ 
We assume from now on that k, p, and s are such that 

90 7^0 

This also implies that 7^ and k_ 7^ 0. 

The function $ is defined in (12.271) up to a muhiphcative constant. In what follows 
we often choose a normalization such that the product of all three branches of $ equals 
one: 

$o$i$2 = 1. (2.30) 

Therefore the function $ with normalization (12.301) is related to the function $ from 
Proposition 12.41 by 

^ ^^^^ 



1/3 ' 
% 

where we recall our assumption that go 7^ 0. Taking into account the explicit expression 
(I2.28P for the equation above, we see that the three branches of the normalized function 
$ behave at infinity as 

'<l>j{z) = -^ + --- , j = l,2, 

z^oo, (2.31) 



where 



LyQZ 

1/3 1/3 

_ % r-< _ % n — 

L^l — ; '-^2 — ; <-^0 



From the equation for $ we can obtain a parametrization of the curve 

Y = {z: |$j(z)| = \^k{z) \ for some < j < A; < 2} (2.32) 
in terms of the function 

J(zy, z) = v^ ^ A{zy + B{z)u + C{z) (2.33) 



where 



A{z) 



3qo - qi{z)q2{z) 

qo 



^^^^ ^ QoqKz) + qUz) - 5qoqi{z)q2iz) + 3gg 34) 

I0 



C{z) 



'^q^qliz) - ql{z)ql{z) + 2ql{z) - Aqoqi{z)q2{z) + ql 



and go, Q'l, Q'2 are the coefficients (I2.28P of the equation for 

Proposition 2.5. The set T given in (I2.32p for the function $ given in (I2.27P satisfies 
V = {z: J{u, z) =u^ + A{z)u^ + B{z)u + C{z) = for some u E [-2, 2]}, (2.35) 
where A, B, and C are given by (I2.34p . 

The proof of Proposition 12.51 is also given in Subsection 13.1.11 
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2.2.3 Structure of F. Definition of cases I and II 



Here we shall use the global structure of F to define the geometrical cases I and IF Since 
the polynomial J(z/, z) in ( ]2.33p is of degree 6 in the variable z, the contour F consists of 
six trajectories zjlu) parametrized by z/ G [—2, 2]. When z/ = 2 we have 



J{2,z) 



U^{z)[P,{z){3s + k) + 9spf 



These trajectories start from the points ai, 6i, a2, 62 and two trajectories start from the 
point 

a = c^ -. 

2>s + k 

Here we assume that a is different from Oj, 6j, j = 1, 2. The case when a coincides with 
one of the aj,bj, j = 1,2 will be treated in the next subsection (see the case HI). We 
denote by 

Tan Tbu Ta2 ) Tfe ' Tan T02 (2.36) 

these trajectories which are then continuously extended as u decreases from 2, see Figure 
12.11 When = — 2 we have 

J(_2, z) = [g0-gl(^)g2(^)]^ ^ ^^^^^ _ ^^y^^ _ ^^y^^ _ ^^y^ ^2.37) 
% 

These six trajectories therefore meet pairwise at the points Pi, (32, Ps, which are the zeros 
of go ~ Q1Q2 (see Figure | 
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bi 







a2 
• 


\/33 
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• 

ai 


• 

bi 


• 

62 



Figure 2.1: Start {u = 2) of the trajectories Figure 2.2: Finish (z/ = —2) of the trajec- 
of F tories of F 



Definition 2.6. We say that the set of points {cti, &i, ^2} belongs to the geometrical 
cases I or II 

A = {ai,bi,a2,b2} elUll (2.38) 
if there exist coordinates B = {k,p, s, c) which are mapped by (12.261) to A such that 
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1. the algebraic function z(z/) defined by the equation J{y^z) = 0, see fl2.33p . has no 
branch points on (—2,2), i.e., that by analytic continuation the trajectories (12.361) 
are defined globally for u G [—2,2] 

7ai('^), lbi{j^), la2{'^), lb2{'^), lai{j^), laii^^), ["2,2]. 

2. When v = —2 we have 

7a,(-2)=7fc,(-2), ^„^(-2)=7,,(-2). 

Now, for y4 G I U II we can define two arcs in C 

=7a,,fe, =7a, U7fe,, J = 1,2, (2.39) 
each connecting aj and bj j = 1,2, and a closed analytic curve 

la • 7ai ' 7^2 • 

Definition 2.7. Given A G I U II, we say that (see Figures 12731 and [2^ 

1. A G I if 7i n72 = 0. 

2. A G II if 7i and 72 have two points of intersection 

7i n72 = {Ci,C2}. 



In the case II we assume that the branch points aj, bj and the intersection points Cj are 
labeled as shown in Figure [27^ That is, if we follow the curve •jj starting at aj {j = 1, 2) 
then we first meet ci and then C2. 





a2 

i 




72 












Figure 2.3: Global trajectories, case I Figure 2.4: Global trajectories, case II 
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2.2.4 Riemann surface for case I. Definition of the global branches for the 
algebraic functions h and $ 



We denote for the case I (see Figure 12. 3p 

■= ■= 7i, J = 1, 2, Ao := Ai U As- 
Definition 2.8. For a set of points A G I the corresponding Riemann surface 



(2.40) 



D^{A) = IHo U fHi U 9^2 (2.41) 
is formed by glueing the sheets of the complex plane, cut along the arcs Aj in ( 12.39^ 

9^, = C\A„ J = 1,2, 

to the sheet 

9^o = C\(AiUA2), 

so that the positive (negative) side of the cut on one sheet is identified with the negative 
(positive) side on the neighboring sheet as in Figure 12. 5[ The curves are oriented from 
aj to bj and the positive side is on the left, the negative side on the right of the oriented 
curve. 





1 1 




' "2 ^^^^^^ 




1 






















An J 






Figure 2.5: Riemann surface for case 1 

The function h in (12.211) is therefore a rational function on the Riemann surface 9^(A) 
in (12.411) with divisor (ll.22l) - (ll.23p . that is, h has simple poles at the branch points a^, 6j, 
simple zeros at the points at infinity on each of the three sheets, and a fourth simple pole 
at a value whose projection is the parameter c. The local definition of the branches of h 
at infinity (11.211) can be extended globally as 

/ioe//(C\(AiUA2)), 

h,eH{C\A^), J = 1,2. 

The function $ in (12.271) is also a rational function on 9^(A) with a double zero at the 
point at infinity on the sheet and a simple pole at the points at infinity on the other 
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two sheets, see f l2.3ip . and normalization (12.301) . The local definition of the branches of 
$ at infinity (12.311) can be extended globally as 

Gi/(C\A,), J = 0,1,2. (2.42) 

The contour F can now be written as a union F = Fo,i U Fo,2 U Fi^2 where 

F,-fc = {zeC \ \^,{z)\ = \^k{z)\}, 0<j<k<2. (2.43) 

In (12.431) it is understood that for z on one of the cuts Ai and A2 we should take limiting 
values. If equahty holds on one side of the cut then equahty holds on the other side as 
well. 

The three branches <l>j, j = 0, 1, 2, also determine a number of regions in the complex 
plane 

fi.-M = {zeC: \<^,iz)\ < |$fc(^)| < |<l>,(^)|}, J,k,i = 0,1,2. (2.44) 
We also define 

n,^k = {zeC: |$,(^)| < |$fe(^)|}, j,fc = 0,l,2. (2.45) 



The contours F^ ^ and the regions ^j^k,e give a partition of the complex plane. 




Figure 2.6: Possible partitions of C by F for case I 



Proposition 2.9. For A we have 

1. There are open sets Uj such that 

Aj c Uj, Uj \ Aj c Qqj. 

2. There exist sets A G I such that (see Figure l2T6ll 

fii,o 7^ or ^2,0 ^ 0- 

Figure 12.61 illustrates possible partitions of the complex plane by the regions and 
the contours Fj^fc. In the left-most figure we have that flifi and f22,o are both empty. In 
the other two figures the set Qi^ is non-empty and it contains part of A2 (as in the middle 
figure) or all of A2 (as in right-most figure) in its interior. 

Now we introduce measures which will describe the limiting behavior of the Hermite- 
Pade approximants. 
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Theorem 2.10. For A in case I we have that the jump of the function Hq on Aq = A1UA2 
produces a positive measure A of total mass 2 

dX{i) = ^ ((/io) + (e) - (/^o)-(O) di, i e Ao. (2.46) 

The measure A consists of two measures of mass 1 with densities (with respect to complex 
line element d^) 

^^■^^^ = 2^K^^^-^°-(^0L/ ^2-^^) 

each supported on Aj, 

Ai on Ai, 



A 

and 



A2 on A2, 



X'(a = ^^M= 7 = 12 

where rrij is analytic on Aj for j = 1,2. 

The proof of Proposition 12.91 and Theorem 12.101 is in Subsection 13.1.21 

2.2.5 Riemann surface for the case II. Definition of the global branches for 
the algebraic functions h and $ 

In the case II we have 71 fl 72 = {01,02}. We assume that the points aj, bj and oj are 
labelled as shown in Figure [231 We use ■jaj,ci to denote the part of the arc •jj between aj 
and Ci and similarly for 7f,^,c2- Then we denote for the case II 

~ ( E-t ■= loX Ao 

A, := 7a,,c, U 76,,C2, j = 1, 2, <^ ' ■ ; (2.48) 

[-£^2 := 7i \ ^1- 

The arcs Ei and E2 form a boundary of a lens- shaped domain G 

dG := El U E2. (2.49) 

Note that the analytic curve 7^ has to pass through the points 01,02 and that it divides 
the domain G into two parts (otherwise it would give a contradiction with the maximum 
principle for harmonic functions). We set 

Ai,2:=7anG. (2.50) 

Finally, we denote 

Aj := A^. U Ai,2, Ao := Ai U A2 = Ai U Ai,2 U A2. (2.51) 
Definition 2.11. For a set of points A G II the corresponding Riemann surface 



m{A) := 9^0 U IHi U 9^2 
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is formed by glueing the sheets of the cut complex plane 

'9^1 :=C\(AiUEi^ 

:=C\(A2UEi) = C\72 

to the sheet ^ 

9=lo:=C\Ao = C\(AiUA2), 

and along Ei the sheets fHi and 9^2 are glued to each other (see Figure [27r|) . 



(2.52) 




Figure 2.7: Riemann surface for case II 

Remark 2.12. 

1. The defined Riemann surface possesses a certain non-symmetry with respect to the 
pairs {ai, hi} and {a2, &2}- We also can use a dual 91 given by 

9^o:=C\Ao, fH2:=C\(A2UE2), 9^i:=C\7i. 

2. Although all three sheets are glued together at the points Ci,C2, it can easily be 
checked that these points are not branch points of 91. 

3. Note that the 9li sheet is a disconnected set. It consists of two components: a 
domain Gi bounded by Ei and Ai_2 

dGi := El U Ai,2, 

and the domain C \ (Gi U Ai). 
The structure of the sheets (12.521) defines the global branches of the functions h and 

$: 



/io,$oe//(C\A^o), 

hi, $1 G H{C \ (Ai U G^)) U H{Gi) 

/l2,$2 Gi^(C\(A2UEi)). 



(2.53) 
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More precisely, in the domains 

C\Ao, C\(AiUGr), C\(A2UEi), 

the branches 

(/lO,$o), (/i2,$2) 

are respectively the result of analytic continuation of ( ll.2ip and ( I2.3ip from point at 
infinity and the branches [hi, $i) in Gi are the result of analytic continuation of {ho, $o) 
through Ai 2 or (/?.2, $2) through Ei. Using the continuity of the global branches of $ 
along 7i, 72, 7a and the maximum principle, we obtain a partition of C by F into domains 
^j,k,e (see f l2.44p for the definition), as is shown in Figure [Z8l 




Figure 2.8: Partition of C by F into ^j,k,e- The case II 

Proposition 2.13. If A e II then the contourV = Fo,iUFo,2UFi,2 defined m ([532]) 
has the following structure: 



ro,i 


:= {z: 




= $1 


:^)i} 


= Ai, 


ro,2 


:= {z: 




= |$2 


:^)i} 


= A2, 


ri,2 


:= {z: 




= |<^'2 


:^)i} 


= El U ^2 U (7, \ Ai,2 



and for the domains ^j,k,i we have (see Figure [2l8]) 

1^^^0,1,2 = 7a UAiU El UE2, 
\dno,2,l = (7a \ Ai,2) U A2 U E2. 

The proof of Proposition 12.131 is in Subsection I3.1.3[ as is the proof of the following 
Theorem 12.141 which introduces the measures. 

Theorem 2.14. For A ^11 we have: 
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1. The jump of Hq on Aq produces a positive measure X of total mass 2 

The measure X consists of two measures Xi and X2 supported on Ai and A2 



A 

with respective densities 



Ai on Ai, 
A2 on A2, 



A'2(0 = , ""^^^^ =, m2 G i/(A2). 

2. The jump of hi on Ei produces a positive measure fii 

- /ii-(O) d^ =: e e i^i 

and fi[ e H{Ei). 

3. There are connections among the total masses of these measures: 

|Ai| + IA2I = 2, |Ai| — = 1. 

Thus the Riemann surface for the case II produces a system of three positive measures 

supp(Ai) = Ai, 
{Ai, A2, /ii}, { supp(A2) = A2, 



supp(/ii) = E 



1) 





+ IA2I 


= 2 


IIA.I 


- l/^ll 


= 1 



If we consider the dual Riemann surface (see Remark I2.12[ item 1), then we arrive at a 
dual system of three positive measures 



(2.55) 





fsupp(A2) 


= A2 


{A2,Ai,/i2}, 1 


1 supp(Ai) 


= Ai 




[supp(/i2) 


= ^2 



and we have 

Ai + A2 = A2 + Ai = A. 



f|A2| 


+ IA1I 


= 2 


\IA2I 


- I/^2| 


= 1 


A. 
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2.3 Geometry of the problem. Case III 

Recall that the geometrical case III is such that the zeros of Pft accumulate on two disjoint 
arcs which do not contain all the branch points. We assume that the branch points are 
numbered such that one arc connects 02 and 62 and such that the other arc connects ai 
and h* with h* 7^ hi, see Figure II. 31 We again start from the algebraic function h in 
fl2.2ip . however we now use its reduced form (ll.24p . Without loss of generality for the set 
A = {ai, &i; a2, ^2} we assume that 

|ai - 6i| > |a2 - 62I, dist(ai, [02, 62]) > dist(6i, [02, 62])- (2.56) 

For a set A with conditions ( 12.56^ we set 

A' = {ai;a2,b2} (2.57) 

and associate with the triple A' the algebraic function f ll.24p : 

^ 01 + 02 + 62 

- S- ^ — h + ^ — = 0. (2.58) 

[z - ai)[z - a2){z - 62) [z - ai){z - a2){z - 62) 

We recall from fll.24p that h in fl2.58p has four branch points at 

/ai + 02 + 62 A ^ 
010262 — I I 

ai, 02, 62, b* = ^ (2.59) 

, , , ^ I ai + a2 + b2 
0102 + ai02 + 0262 — 3 



We do not need to use the coordinates fl2.23p now since we already have the explicit ex- 
pressions in terms of the input data A' for the coefficients of the equation (12.580 . However, 
in order to avoid too cumbersome expressions we set, without loss of generality, 

ai := —1, 62 := 0, 02 := a, (2.60) 

and a 7^ —1, a 7^ 0. 

Proposition 2.15. The exponential function of the Abelian integral 



of the function h in (12.580 is, up to a multiplicative constant, an algebraic function satis- 
fying the equation 

$^ + gi(^)<l>^ + g2(2)$ + go = 0, (2.61) 

where deg qj < j for j = 0, 1, 2 and the qj are rational functions of ai, 02, 62- If we take 
into account (I2.60p then we have 



2 , 5a ^ , a(a^ + 4a + 1) 



qi{z) = z(a-l)(a' + y + 1) + 

/27z^ 9 1 
q2{z) = -K\^-^--{a-l)z--{a' + Wa + l)] (2.62) 

go = n'^ 



where k = a {a + 1) /4. 
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The proof of Proposition 12.151 is in Subsection I3.1.4[ 

Now substitute the qo, qi, q2 from fl2.62p into the expressions (12.341) for the coefficients 
A,B,C (see Proposition 12.51) of 

J(zy, z) = u^ + A{zy + B{z)v + C{z). 

Then as in Proposition 12.51 we obtain a parametrization of the contour F 

T := {z : \^j{z)\ = |$fc(z)| for some < j < A; < 2} 

= {z : J(z/,z) = for some z/ G [-2,2]} (2.63) 

as the union of six trajectories. When u = 2 the three trajectories jai, 7a2) 7^2 start from 
the points Oi, 02, 62 and the other three trajectories 7^^^ {j = 1, 2, 3) start from the point 
b* (see Figure [2^ . When u = —2 these trajectories meet pairwise at the points Pi, P2, P2 
given by (12.371) (as in Figure [2^ . 





0-2 
I 




A \ 















Figure 2.9: Start [v = 2) of the trajectories of F for Case III 

Definition 2.16. We call the triple A' in (12.571) acceptable for case III if 

1. The algebraic function 2;(z/) defined by (I2.63P has no branch points on (—2,2); 

2. When 1/ = —2 we have 

1a,{-2) = lb,{-2), 7a,(-2) = 7j'V2) for some j G {1,2,3}. 

If part 2 of the definition is satisfied we assume without loss of generality that the 
trajectories starting from b* are numbered such that 

7a.(-2)=7iP(-2). 

It follows from Definition 12.161 that the trajectories 7a2, 762, 7^^, 7^!'' for an acceptable 
triple A' are defined globally for u G [—2,2] and we can define the arcs joining the branch 
points (I2.59P as 

At:=7aiU7jP, A2:=7a2U7f,2, Aq := U A2. (2.64) 

Then we define for an acceptable triple A' the Riemann surface as in subsection 12.2.41 (see 
(1231])) 

m*{A') = 9^* u in* u (2.65) 
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with three sheets 

1H* = C\(A*U As), mi=C\Al, % = C\A2, 

and the global branches of the algebraic functions h and $, defined by ( I2.58P and (]2.6ip . 

are 

ho, $0 e H{C \ {Al U A2)), hu $1 G H{C \ A*), h2, $2 e H{C \ A2). (2.66) 

Possible partitions of C by F into the domains ^j^k,e are shown in Figure 12.101 We 
have 

Proposition 2.17. 

1. There exists an open set U2 such that 

A2CU2, ?72\A2C 1^0,2 

and there exists an open set Ui such that 

At\{6*}ct/i, t/i\At cfio,i. 

2. For any acceptable triple A' we have 

1^1,0^0 and h* e BVLq^i. 

In addition we have that VLi q is connected. 

Moreover the trajectories 'j^p {j = 1,2,3) start from b* at an angle 2^/3 and split a 
neighborhood of b* into two domains belonging to Qq i (the boundaries of the domain 
contain Al) and one domain belonging to Qi^ (the boundary contains %J and %J )■ See 
Figure^^ 

Now we introduce the measures. 

Theorem 2.18. For an acceptable triple A' we have that the jump of the branch ho on 
A = AJ U A2 produces a positive measure X on '-f of total mass 2. The measure consists 
of two measures of mass one: 

f^'i(0 = 2ii[(M+(0-(M-(0], ^eAi, 

A'(0 = 

U^(0 = 2^[(M+(0-(M-(0], eGA2, 

and 

V e - ai V(^-«2)(e-&2) 
with mi analytic in a neighborhood of Al and m2 analytic in a neighborhood 0/ A2. 

The proofs of Proposition 12.171 and Theorem 12.181 are in Subsection I3.1.4[ 
We can now define the geometrical Case IIF 
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^0,1, 



^0,1,2 



Figure 2.10: Partition of C by T. The case III (acceptable triple) 



Definition 2.19. We say that a set of points A := {ai, bi, 02, 62} satisfying (I2.56p belongs 
to the geometrical Case 111 

A e 111 

if 

1. The triple A' = {ai, 02, 62} is acceptable. 

2. bi e Qifl. 

For A G 111 we assign to A the Riemann surface IH* and the algebraic functions h and 
$ corresponding to its triple A', i.e., we use fl2.58p . (12.611) . and (12.651) . 

2.4 Geometry of the problem. Common branch point: cases IV 
and V 

In this section we present a complete classification of the geometry of the Hermite-Pade 
approximation problem for the two functions (11.31) with a common branch point, i.e., 

fj eAiC\{a,,b}), ai^a2. 

i.e., 

A = {ai,b] a2,b}. 

Again, as in the case 111 we associate with the triple 

A' = {ai, a2, b} 



(2.67) 



the algebraic function h defined in (I2.58p . and the Abelian exponential $ given by (I2.6ip . 
There are three possibilities for the position of the points A: 



1. The triple A' in (I2.67P is acceptable (see Definition 12. 16p and b G fii^o- Then A can 
be realized as a limiting situation of case 111 and we say 

A e 111, 

(see top figure of Figure [2711!) . 

2. The triple A' is acceptable but b ^ (see Figure I^TTT] in the middle). 

3. The triple A' is not acceptable (see bottom figure of Figure [2. lip . 



37 



Ill 














b*/ 






ai 


^0,1 









V 






V 

Y' ^1,0 


ai 


^0,1 







IV 




c\ 


-«2 1 

( V* 




^0,1 




\ A ^0,2,1 

^6 J 



Figure 2.11: The subclasses III, V and IV for a common branch point 

Definition 2.20. We call the subclasses 2 and 3 above the cases V and IV respectively. 

As in the previous subsection (see (12.601) ) we set, without loss of generality, 

a\ = —1, 6 = 0, 02 := a, \a\ < 1. 

In order to give the explicit description of these three subclasses we define three regions 
Dj of the disc D = {a : \a\ < 1} as follows 

aeDi A = {-1, 0; a, 0} G III, 

aeD2 AeV, (2.68) 
a e A e IV. 

We have 
Theorem 2.21. 
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1. Two branches of the algebraic curve 

+ (4 - 4z/)a^ + (22 - 8u)a'^ + (4 - 4z/)a + 1 = 
lying in D for v G [—2, 2] forra the boundary dDi of the domain Di. 

2. The algebraic curve 

P{a, v) := 16a^2 ^ gg^n _^ (335 _ i08zy)a^° + (800 - 540z/)a^ + (2169 - 1404z/)a^ 
+ (4932 - 2376z/)a^ + (6630 - 2808z/)a^ + (4932 - 2376i/)a^ 
+ (2169 - 1404z/)a^ + (800 - hAQu)a^ + (336 - lQSu)a^ + 96a + 16 = 0, 

which for u = 2 is factorized as 

(a-l)^(2a + l)^(a + 2)^ = 0, 

has six branches in D for v G [—2,2]. These six branches start (for u = 2) from 
the points 1 (two branches) and —1/2 (four branches) and give the outer boundary 
dD'^^ of the region D2 ■ The boundary dD2 therefore consists of 

dD2 = 9Di U9D°"*. 

3. The region is the open set bounded by dD"^^ and {\a\ = 1}. 
The proof of Theorem 12.211 is in section 13.1.51 
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Figure 2.12: The regions D2, -D3 and their boundaries 
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Now we have to assign the Riemann surface to these subclasses and we have to define 
the global branches of the algebraic functions h and The subclass III is already consid- 
ered (see subsection I2.3p . For the subclass V we again use (as in III) all the attributes of 
the acceptable triple A', i.e., the Riemann surface fH*(A') in ( 12.65^ and the global selection 
of the branches (12.661) . For the subclass V we can also make item 2 of Proposition 12.171 
more precise. 

Proposition 2.22. For A we have (see Figure d. 11\ second picture) that b G ^o,2,i- 
The proof of Proposition 12.221 is in section I3.1.5[ 

We denote the connected component of flo^2,i containing b by Gi and its boundary by 
El. Also we set (see (EJl])) We set 

Ai,2 := AsHGi. 

For the subclass V we use the measures Ai, A2 introduced in Theorem 12.181 We set 

Ai := Ai, Ai,2 := A2 , A2 := A2, 

Al,2 

where Ai and A2 at the right hand sides are from (12.661) . In addition to these measures 
we shall use the balayage of the measure Ai,2 on the contour Ei. We denote this measure 
by Hi. 

^w = -^Ai,2^ on El. (2.69) 

The last case IV requires special treatment. We start with a characterization of the 
contour T = {z : \^j{z)\ = \^k{z)\,j 7^ k,j,k = 0,1,2} for a non- acceptable triple A'. 
Consider the alternative to Definition 12.161 of an acceptable triple (see Figure 12.131) 

Definition 2.23. 

A. If the algebraic function z(z/) has a branch point on (—2, 2) (see (I2.63P ) then the triple 

A' is called critical. 

B. 'The triple A' is called strictly non-acceptable when 

1. The first condition in Definition 12. 161 still holds. 

2. For z/ = 2 we have (maybe after renumbering of trajectories 7^^'') 

7.,(-2) = 7^P(-2), 7a.(-2) = i!;\-2), ^,{-2) = ^S\-2). 

From the definition we derive (see Figure I2.13p 
Proposition 2.24. 

A. For a critical triple one of the trajectories •ja, « £ {^2, b} meets the trajectory 7^, or 

(3) . 

3z/^ G [-2, 2), 3a G {02, h}, 3j G {2, 3} : 7„(z/^) = -f^^\iyj. 
For a critical triple the meeting point of the trajectories coincides with b* 
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B. For a strictly non-acceptable triple the arcs 

7i = 

intersect at one point c 



I I (1) I I (2) I , (3) 



(2.70) 



{c} = 7i n 72 n 73 7^ 0. 
The proof of this proposition is also in Subsection I3.1.5[ 
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Figure 2.13: Non-acceptable triples: left is A-critical, right is i?-strictly non-acceptable 

For strictly non-acceptable triples we now assign a Riemann surface and define the 
global branches for the algebraic functions h and $. We have that each arc from (12.701) 
is divided into two parts by the point c: 



I I (1) I I (2) 11 (3) 

7l = 7ai ,c U 7^'. , 72 = 7a2,c U 7^'. , 73 = 76,c U 7^'. . 



(2.71) 



We take three sheets of the complex plane cut as indicated (see Figure [2m|) : 

= C\(7ai,cU7a2,cU7fe,c) 

= C\7i 

% = C\(7lilU7a2,cU76,c), 

and then we glue them together so that they form a Riemann surface of genus zero: 



2.72) 
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which we assign to the strictly non- acceptable triple A'. The algebraic functions h and $ 
for the set of parameters A' = {ai;a2,b} are single valued rational functions on fH*(y4'). 
The structure of the three sheets (I2.72p defines the global branches of h and $: 



ho,% e i7(C\ (7a,,eU7a2,cU7fe,e)), 

G i/(C\7i), 

/12,<1>2 G i/(C\(7SU7a2,cU76,c)). 



(2.73) 



We denote 



and 



We also denote 



Ai 



Ai := 7ai,, 
Ai U Ai2, 



A2 := 7 



^1,2 



7c,6, 



A2:=A2UAi,2, Ao:=AiUA2. 



E2 := 7i \ Ai, El := 72 \ A2, Eq := 73 \ Ai,2. 



(2.74) 
(2.75) 
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• 




ai 





Figure 2.14: Sheet structure of the Riemann surface for a strictly non- acceptable triple 
(case IV) 



Proposition 2.25. The contour T (see fl2.32l) - fl2.33p ) has the following structure: 



ro,i := {z:\<^o{z)\ = \^i{z)\} = ^u 
ro,2 := {z:\%{z)\ = \^2{z)\} = A2, 

3 

ri,2 := {^:|<l'i(^)| = |$2(^)|} = |j7Jv = ^oUEiUE2, 



(2.76) 



and for the domains ^Ij^m we have (see Figure d,. Ill third picture) 

^0,2,1 = G U G, 

where the components G and G are defined by their boundaries 

dG = EiU E2, dG = EoU E2, 

and 

^0,1,2 = C \ ^0,2,1- 

Now we introduce the measures. 
Theorem 2.26. For A E TV we have 

1. The jump of Hq on Aq produces a positive measure A of total mass 2 



-^(/io+(0 - /^o-(O) di =: dA(0, e G Ao. 



(2.77) 
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The measure A consists of two measures Ai and A2, with |Ai| + IA2I = 2, supported 
on Ai and A2 respectively 

Ai on Ai, 
^ A2 on A2, 



A := 



Ai(0 = ^iilL, miGi/(Ai), 

- «i 

A2(0 = , ""'^^^ =, m2G/f(A2). 
^. r/ie jump of h2 on E2 produces a positive measure fi2 

and 

3. There are connections among the total masses of these measures: 

|Ai| + IA2I = 2, IA2I — |/i2| = 1- 

We also denote 

Ai,2:=A|^^^^, Ai:=A|^^, A2:=A|^^. (2.78) 
The proof of Proposition 12.251 and Theorem 12.261 is given in Subsection I3.1.5[ 

Remark 2.27. [on critical cases] There are several cases of location of the branch 
points A = {ai, b; 02, b} which are the hmits of the cases described above. For example, 
the case (considered by Kalyagin in [H]) when 

Ao:={-l,0;l,0} 

is the limit of the class I 

A,,, := {-1, -e; l,e}, as e 0, 

and at the same time it is the limit of the class III 

^o,e = 0; 1, e}, as e — > 0. 

Also, the critical triple is both the limit of the strictly non-acceptable triple (case IV) and 
of the case V. We call these cases the critical cases. For the critical cases the assignment of 
the Riemann surface and the definition of the global branches for the algebraic functions 
h and $ is carried out by passing to the limit in the non-critical cases described above. 
Even though one of the main ingredients of our analysis, namely the geometry of the 
problem, is clear for the critical cases, we do not consider the Hermite-Pade asymptotics 
for these cases in this paper. (In what follows we exclude the critical triples from the 
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class IV.) The point is that the rigorous proof of the asymptotics for these cases requires 
a non-standard local Riemann-Hilbert analysis. 

Certain critical cases of 2 x 2 matrix- valued Riemann-Hilbertproblems have been stud- 
ied recently and it was found that the local Riemann-Hilbert analysis can be constructed 
in an explicit way with the isomonodromy characterization of certain Painleve transcen- 
dents, see [IHl EH [22l [23l Ell ESI 133] • Similar constructions can be expected to apply to 
some of the critical cases related to this work, but there are other critical cases as well 
which are special for 3x3 matrix-valued Riemann-Hilbertproblems and which do not 
occur for 2x2 problems. 



2.5 Weak asymptotics, convergence and vector potential prob- 
lems 

In this section we formulate corollaries from the strong asymptotics of the Hermite-Pade 
polynomials regarding the weak asymptotics of the poles of the Hermite-Pade approxi- 
mants and their convergence. The weak results follow from the strong asymptotic formulas 
(Theorems 12.351 and 12.361) but they are stated here since to present these results we only 
need the basic functions which we have defined during the description of the geometry of 
the problem. The results will be different for each geometrical case. 

2.5.1 Weak convergence 

We recall that we consider the Hermite-Pade approximants for two functions (12.191) . see 
Definition 12.11 and convention ( 12. 16^ : 

/, e ^(a„ a,; 6„ /?,;(],), j = 1,2, (2.79) 

where the domains Qj for the analytic continuation of the weight functions wqj, defining 
fj (see (12. 5p and (12.21) ). depend on the location of the branch points 



A = {ai,6i;a2,&2}- 

We assume flj is such that: 



1) njDAj, J = 1,2, for A el 

2) D Aj UG, j = 1, 2, for A e II, (2.80) 

3) l^' ^ for A E III, 

4) nj D Aj UG, j = 1, 2, for A e IV. 

5 <^ ^ ' for AG V 

1 1^2 D As U Gi, 



We recall that the definition of the classes I, II, HI, IV, and V is given in Definitions 12. 7[ 
12.191 12.201 and the corresponding arcs Ai, A*, A2 and sets Gi, G are defined in (I2.40p . 
(ICT]) . dlMD, Proposition EJl (I^TD . and ([2Tri) . 

The weak limit of the counting measures z/p^, which have equal mass l/{2n) at the 
poles of the Hermite-Pade approximants (see (11.41) ). has a universal character. 
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Theorem 2.28. Suppose that A belongs to one of the cases I, II, III, IV, or V. Then the 
poles of the Hermite-Pade approximants (11.21) for the functions in (12.791) - fl2.80p have a 
weak limit 

t/p„ A A/2, n ^ oo, (2.81) 

where the limiting measure A is defined in Theorems \2.1(\ \2.14\ \2.18[ \2.26\, depending on 
the geometrical class. 



Now we state a result about the finite zeros of the functions of the second kind (11.21) 

Rl^^ := f,P„- Q^^\ J = 1,2. 

These zeros represent the extra interpolation points. We use the notation u for the 
counting measures with equal mass l/n at the finite zeros of Rn\ 

Theorem 2.29. Consider Rn^ for the functions l[2. 79^ where f2j is as in l{2.80\) . 

1. When A G 1 U 111 there are no extra interpolation points for n large enough. 

2. When A G 11 U IV we have 

J = 1,2. (2.82) 



tin 



/ij, 



(2) 

3. When A G V the function Rn has no finite zeros for n large enough and 

z/ (1) A 



In 2 and 3 the limiting measures /ii, /i2 are defined in Theorems \2. 141 \2.2(k and in \2. 55]) . 
and Ii2. 69\) . depending on the geometrical class. 

The next theorem describes the nth root asymptotics of the error term 



fj - ^ 



U) 



71. 



IT. 



P ' 



J = 1,2, 



and the convergence of the Hermite-Pade approximants (11.21) . 



Theorem 2.30. Consider the Hermite-Pade approximants for the functions ^2. 79 ) with 
Qj as in ( f^.gOj) . 



1. When A & l we have 



(j-)|l/n 



uniformly on compact subsets ofC\Aj, j = 1, 2, and therefore for the Hermite-Pade 
approximants we have 



TT, 



(i) 



on fii 



oo, on C \ (Aj U f^oj), i = 1, 2, 

with geometric rate. Furthermore (see Figure \27^) there exists an A such that 
^1,0 ^ or 1^2,0 7^ 0- 
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2. When A G III we have 



If. 

\h "n I 



uniformly on compact subsets of C\ when j = 1 and C \ A2 when j = 2, and 
therefore, with the convention 112. 56]) . we have 



7rW-./i on C\ (At ur^i.o) 

kn ^ - /il ^ 00 on nifi, 



IT, 



(2) 



/2 on C\A2 



Furthermore (see Figure \2. 1 (A) Vti^ ^ 0. 
3. When A&V we have (uniformly on compact subsets of the indicated sets) 

on (C\ A*) \Gi 

on Gi \ A2 



/■ 


$0 


< 






$0 




$2 



r2)|i/n 



on C \ A2 



1/2 -vri^ 
anc? therefore (see Figure 12J1\) 

onc\(Atuni;^) 



where Qi q 7^ 0, and 



- /il ^ 00 on fii,o, 



42)^/2 on C\A2. 



^. When A ^ IV we have (uniformly on compact subsets of the indicated sets) 

on (C\ Ai) \G 
on GVAi^a 



/l - 





$0 


< 


$1 




$0 




$2 



I/2 -71", 

and therefore we have 



(2) 1 1/1 



$2 



on C \ (A2 U E2 



f, on C\Aj, j = l,2. 



5. When A G II we have (uniformly on compact subsets of the indicated sets) 



/l - 



$0 
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on C \ (Ai U El 



1/2- ^''^ 





$0 


< 






$0 




$1 



on (C \ A2) \ G 
on G\Ai,2 



and therefore we have 



IT, 



f, on C\A„ J = 1,2. 



We recall that the definition of the branches of the algebraic function {$0; '^'i; '^'2} and 
the domains of divergence ^jfi, j = 1,2, depending on the geometrical case is given in 

dllS]), dMZD-dSSHD, dm, dSl), dMlD, and 



2.5.2 Vector equilibrium problem 

Concluding this section we state a universal vector equilibrium problem for the loga- 
rithmic potentials of the measures A = Ai + A2, /ii, /i2- These measures were introduced 
in Theorems 12.101 12.141 12.181 and 12.261 for the description of the weak limits of the poles 
and of the extra interpolation points (see Theorems 12.281 and 12.29^ . 

Theorem 2.31. Suppose that 



A G I U III U II U V U IV, 



then 



I. a) There exist piecewise analytic arcs Ai, A2 which make the functions /i,/2 in 
^2. 79{ ) holomorphic 

/, Gif(C\A,), J = 1,2, 
and a piecewise analytic contour E containing the common part of Ai and A2 

Ai,2:=AinA2, (Ai,2 = ^ ^ = 0). 

h) There exists a triple of measures (Ai, X2, fJ'i) with supports (we use the notation 
S{fi) = supp(/i); 

Six,) C Ai, SCX2) C A2 := A2 \ Ai,2, 5(/ii) C E, 

and with the relations on their total mass 

|Ai| + IA2I = 2 
|Ai| - = 1. 

c) This triple of measures possesses the following equilibrium relations with some 
constants ni and /t2." 



2V^^ + V^^ - 



= Ki, on ^(Ai), 

> Ki, on Ai, 

= K2, on S{X2), 

> k2, on A2, 

: K2 — Ki, on S{fj.i] 
> 1^2 — Hi, on E. 



(2.83) 
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f dUi 


dUi 




dn^ 


dU2 


_ dU2 


dn^ 


9n_ 


dUs 


_ dUs 


. dn+ 


dri- 



d) The supports of the measures possess the following symmetry relations 

on ^(Ai), 

on ^(Aa), (2.84) 
on S{fii), 

where d/duj. denotes the normal derivatives on the respective contours. 

II. There is also a dual problem regarding the triple (A2,Ai,/i2) which can be obtained 

from the problem I.a)-I.d) by interchanging the indices 1 and 2. 

III. The equilibrium measures (Ai,A2,/ii) and (A2,Ai,/i2) are related as follows 

A := Ai + A2 = A2 + Ai, 5(A) C Ao:= A1UA2, 

/x:=/ii + /i2, ^(yUi) U S'(/i2) = -E, 
yf^ =y^i,2 ^ Ai,2:=A 



IV. The measure A/2 zs the weak limit h2.81\) of the poles of the Hermite-Pade approxi- 
mants of the functions \2. 19^ , and the measures fii and ^2 are the weak limits h2.82\) 
of the extra interpolation points. 



The proof of Theorem 12.311 is given in Section 13.21 

We now describe the basic notions of the equihbrium problems I and II of Theorem 12. 301 



Ai,A2,^ 
Ai, A2, /ii 



A2,Ai,E 

A2, Ai, /i2 



in conformity with the geometrical cases. 
Remark 2.32. 

1. In the case I we have Ai^2 = and therefore 

'a,- = A„ j = 1,2, 
^ = 



(2.85) 



Aj = Aj, yUj = 0, S'(Aj) = Aj, j = l,2, 

and in the systems of the equilibrium and symmetry relations (I2.83l) - (l2.84p only the 
first two relations are needed. 

2. In the case III, we set 

Ai = A* U 51, 51 C fii,o (2.86) 

where 51 is an arbitrary rectifiable arc in i7i,o,2 joining the points h* and hi (see 
( I2.45P and Figure [2.10p . For this case we have the same setting (12.850 as for the 
case I with the minor difference that S{\i) = A* C Ai. 
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3. We note that for the cases I and III the equihbrium relations (12.831) reduce to the 
equihbrium relations for an Angelesco system (11.121) . (ll.lSp in the complex plane. 

4. In the case V we set (see Proposition 12.241) 
Ai := Al U 61 U Ai,2, A, := A, \ Ai,2, j = 1, 2, 

where 5* is an arbitrary rectifiable arc in f2i,o,2 and joining the corresponding end 
points of the arcs A* and Ai 2- Here we have 

E = El ^ 0. 

If we consider the dual problem (see II of Theorem I2.3ip for the triple 

A2,Ai,E 

A2, Al, /i2 

then the solution of the Angelesco equilibrium problem for the case III 

'2V^^ + V~^^ = K2, onA2, 

_^ 2V'^^ > ki, on Ai(= Ki on AJ) 

provides the solution of the dual equilibrium problem with /i2 = 0, because (see 
Figure I2.1H second picture) 

V^^ -ki = V^^ - K2, on E. 
In order to obtain the measure [i\ we have to take the balayage of the measure 

Al 2 := A2 

Ai,2 

on the contour E or consider the problem (I2.83p - (l2.84p for the triple 

Ai,A2,E, 
Al, A2, 

5. In the cases IV and II the contours for the equilibrium problems are defined in 
^HM-^^ and ([23H])-(E3ID- In these cases (see Theorem [211 and OS]) there 
are no degeneracies of the components of the equilibrium problem (I2.83P and the 
dual problem. 

2.6 Szego functions for the Hermite-Pade polynomials 

The Szego function is an important ingredient of the strong asymptotic formulas for 
orthogonal polynomials on an interval of the real line and on the unit circle. In this section 
we introduce a generalization of the Szego function which we need for the presentation 
of the strong asymptotics of the Hermite-Pade polynomials. We shall define the Szego 
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function as a solution of a certain boundary value problem on the three-sheeted Riemann 
surface 

2 

j=0 

Before we state the boundary value problem we need some preparation. We fix the cuts 
Sj^k (Jordan arcs in the complex plane or a union of Jordan arcs) where the sheets Dij and 
dlk are glued together 

fH,=:C\(5,,,U5,,fc), i^{k,j}, k^j, i,k, J =0,1,2. 

The two sides of the cuts (5^ j U Si^kY'^'^ form a boundary of the sheet (Jordan curves 
of the Riemann surface 9^) 

dm^ := dm^j U dD^e,k, dD^e,k := % H = 5^ J U 6^ . 

Thus 

in =: U dm, dm := (9fHo,i U (9fHo,2 U 91Ki,2. 

We fix the orientation of the union of the Jordan curves dm such that 

99^5-^ c M„ dm^^r^^ c Mfc, J > k, J, k = o,i, 2. 

We also fix an orientation of the Jordan arcs 6j^k- Our Riemann surfaces have four branch 
points at {ai, bi, 02, 62}, with a possible replacement of bi by b* for the cases III, IV and 
V. We have 

1. {aj, bj} C Soj, J = 1, 2 for A e I U II, 

2. {ai, b*} C 5o,i, {^2, C 5o,2 for A G III U V, 

3. ai G 5o,i, 6* G 5i,2, {02, &} C ^0,2 for A G IV. 
Let ujj{z) be the branch of the function 

uj'^{z) = {z — aj){z — bj), ujj{z) = z + ■ ■ ■ , z ^ 00, j = 1,2. (2.87) 

For the cases III, IV and V we replace bi by b* in fl2.87p . From the statement of the 
problem (see Section [2?T|) and from the definitions of the geometrical cases (Sections 12. 2f - 
12. 5p it follows that the weight functions wi and W2 are defined on the arcs (5o,i and 5o,2- 
We set 

Wj := iWjUj-, on 6oj, j = 1,2. (2.88) 

We assume that the orientation of the arcs 6kj is fixed so that (see, for example. Figure I^TTI) 
each Wj has an analytic continuation on 5i 2 (when 61^2 7^ 0)- Then we define 

{Wi, on 5o,i, 

W2, on 5o,2, (2.89) 
W2/W1, on (5i,2. 
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We note that, depending on the position of the points A, we have 

A e I U III U V ^ (5i,2 = 0, 
A G II U IV ^ (5i,2 ^ 0. 

Finally we duplicate the weight function w on both sides of 5^^^ i.e., we define w on dDl 

w+ = w_ on 6kj =^ w on d^i. (2.90) 

Now we formulate the boundary value problem on 9^ whose solution defines the desired 
generalization of the Szego function. We are looking for a piecewise holomorphic function 
JF on IK such that 

'1. T e H{9\\dm), 

2a. 3J-± G C{dD\ \ {ai, a^, h}) : -F+ = T-W on 8^ 

2b. \J^\z)w{z)\ =0(1), z^A ^ ■ ^ 

^3. ^(ooW)^(oo«)^(oo(2)) = 1. 

Remark 2.33. The maximum principle implies that if from fl2.9ip exists, then 

a) For every z G C 

J^(z(°))J^(^(i))J^(^(2)) ^ ^ (2.92) 

b) JF is unique. 

Here z^^^ := 7r^^(z), where vr : — > C is the projection of the sheets of 91 onto C. 

Due to (12.901) the index of the boundary value problem (I2.9ip is equal to zero. This 
gives the existence of the solution of (I2.9ip . In order to write an expression for the solution 
of (I2.9ip we use the meromorphic (Cauchy) differential on 91 

dM^{z), zg91,^g91, 

which has simple poles at the points ^ G 91 and at the points ^ G 91 where ^ has the same 
projection on C as ^, i.e., 7r(0 = vr(f), ^ 7^ f : 

and the residues at the poles are 

resM^ = 2, resM^ = -1, 7r(0 = 7r(0, ^ 7^ I 

Then Cauchy's theorem on 91 and (I2.92p give 

^ [ log*(;2)ciM5(z) = 21og.F(0-log.F(fW)-log.F(f(')) = 31og.F(0. 
Summarising we have proved 

Theorem 2.34. There exists a unique solution of the boundary value problem Ii2.91\) . 
which is given by 

J^iO = exp /^log*^'^) ^^«(^)) ' e e 91 \ 991. 

Later on we shall use the notation for the branches of the multivalued function JF 

Fe{z) := T{z^'^), £ = 0,1,2, (2.93) 
where z eC \ {Sej U Sg^k), ^ ^ {k,j}, k ^ j, and i, k,j = 0, 1,2. 
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2.7 Strong asymptotics 



We can now state the main results of this paper on the strong asymptotics of the Her- 
mite-Pade approximants. The formulas for strong asymptotics contain the functions $ : = 
{$0, $1, $2} from (|T:28|) and the Szego functions T := {Fq, Fi, F2} from (ICTjl and (IZMll . 
These functions are defined by the Riemann surface D\ of the corresponding geometrical 
cases (defined in Sections l2.21 - |2.5p . 

The strong asymptotics of the denominators P„ of the Hermite-Pade approximants 
have a universal character. 

Theorem 2.35. Suppose that 

A = {ai, hi] a2, 62} e I U II U III U IV U V. 

Then the denominators of the Hermite-Pade approximants (11. 2p for the functions fj, with 
branch points {oi, 61, 02, &2} (see (I2.79P - (12.801) ) have the following asymptotic formulas as 
n 00, which hold uniformly on compact subsets of the indicated sets: 



Pniz) 



i"o(oo) 



{Co%{z)) "(l + 0(l/n)), ^gC\(5o,iU5o,2), 



and 

Pn{z) 



Fo{cx>) 



(Co$o(^))' 



+ 



Fojoo) 



(Co<fo(^))^ 



(1 + 0(1)), 
z G (5o,i U 5o,2) \ A. 



We recall that we have 





= A„ 


j 


= 1,2, 


for A el, 


^ ^0,1 


= Al, 




= A2, 


for A G III U V 


^1 


= Ai, 




= A2, 


for A G II, 


Ml 


= Ai, 


So,2 


= A2, 


for A G IV. 



(2.94) 



The strong asymptotics for the functions of the second kind Rn \ j 
forms, depending on the geometrical case. 



1,2, have different 



Theorem 2.36. For the asymptotics, as n —>■ 00, of the functions of the second kind Rn\ 
given by (12.70 . of the Hermite-Pade approximants for the functions fj with branch points 
{01,61,02,62} in (12T9|) - (I2:80D . we have 



1. When A elVMW we have 

Fo(oo) 



ujj{z)Fj{z[ 



.(Co$,(;.))-"(l + 0(l/n)), 



i"o(oo) 



iUj{z)Fj{z 

for j = 1, 2. Here the 5oj ore as in (12.941) 



(Co«f,(.))~"(l + 0(l/n)) 



z G C\5o,i, 

z G 5o,i \ A, 
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2. When A^N we have 

UJ2[Z)l'2[Z) 

^ (^°^^-("))""(^ + 0(V^))) ^ > .eAAA 

and 

= f- ?[p] A CoM^))-") (1 + 0(1/")), ^eA;\A 

zedGi\ A. 

3. When AelV we have, with 5o,2 := U 

z e ^2 \ {&*}, 
anc? wz^/i (j such that dG = EiU E2 

= (- Ji|;y(^o..(.))-"(l.0(l/n)))^^, zeA.,2, 
^ ^oM f I I ](l + 0(l/n)) 

"^^ Q V iOl{z)F,{z)^U^) W2{Z)UJ2{Z)F2{Z)^U^)J^ ^ 

zeE,\{b*}, 
^ Fo(oo) l + 0(l/n) 
" ^ ^ Q Mz)F,{z)<!>-,{z)U 

Fo(oo) l + 0(l/n) , ^ p \ r.*; 

53 



4- When A E II we have the same asymptotic formulas as in 3) (the case A E IV) if 
we interchange the indices 1 and 2 and replace {h*} by 0. 

The proof of the strong asymptotics formulas is given in Section HI It follows from a 
steepest descent analysis of the Riemann-Hilbertproblem (12.131) . 

3 Proof of the geometrical results and equilibrium 
properties 

3.1 Proof of the geometrical results 

Here we present the proofs of the propositions and the theorems stated in subsections 
I2.2H2.4[ We start with some results of general character. 



3.1.1 Proof of Propositions [2l3l EH and [2I5] 

Proof. 

1. The parametrization of the algebraic curve fl2.2ip of order 3 and genus 0, given in 
Proposition 12.31 can be checked directly. Indeed, substituting the polynomial coef- 
ficients (12.241) and (I2.25P into the equation (I2.2ip for the function h and computing 
the discriminant 

_ p| - n^Pl 

we obtain that the polynomial 

P| - n4p2 = V = [kPl + 3Pips + S^f 

is a complete square and therefore its zeros are not square root branch points of 
the algebraic function h (the Puiseux series at these points have integer exponents 
or cubic roots in some degenerate cases). Thus (in the non-degenerate cases) the 
only branch points of the function h are the four zeros of the polynomial 114. The 
Riemann-Hurwitz formula (for a function of the third order with four branch points) 
implies that h has genus 0. 

2. To find the coefficients (I2.28P for the equation of the algebraic function $ in Propo- 
sition [2]1] we proceed as follows. We make a linear change of variable z ^ Pi := z — c 
and substitute the series 

CO 

k=l ^1 

into the equation (I2.2ip to find two sets of coefficients {c[,^^} and {c^^^}. Thus, we 
obtain the coefficients of the expansion at Pi = oo for the three branches of h 

°° Jl) °° ^(2) °o ^(1) , J2) 



k-- 



_j pk ' / V pk ' / J pk 

=1 1 fe=l 1 k=l 1 
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Then, see f l2.27l) . we integrate these three series, take the exponential function of 
the results and again expand it into a power series around Pi = oo 

oo ,(1) oo ,(2) oo ,(0) 

fc=0 1 k=0 1 k=2 1 

Next, we find the constants m^, £ = 0, 1, 2, such that the power series 

2 

is a polynomial. In this way we obtain the series around Pi = oo for the branches 
of the function $ 

= TTii^e. (3.1) 

Finally, we use the Vieta relations to obtain the coefficients (12.281) from the coeffi- 
cients of the series (13. ip . 

3. To obtain the algebraic parametrization of the curve F in Proposition 12.51 we note 
that the coefficients of the polynomial J(z/) 

.\ = ^ ^iM _ ^ (^M : ^iM _ ,^ : £iM _ 

(3.2) 

are symmetric functions with respect to $0)'^i;'^2- Representing these symmetric 
functions by means of the basic symmetric functions defined by the coefficients of 
the algebraic equation for $ we arrive at (I2.33p -( l2.34p . 

□ 

3.1.2 Proof of the geometric results for the case I 

Next we prove the results concerning the geometrical case F 

Proof of Theorem \2.1(A and Proposition \2.9\ . First we notice that (12.460 defines a charge 
of total mass 2. Indeed, 

/ dX = ^[ (V(0-/io-(0)rfe=7^ / MOC 
Jao 27r2 27rz Jg^ 

and by Cauchy's theorem and (ll.2ip we have 

/ dX = 2. 

Analogously we have 

f d\^ = l, J = 1,2. 
Ja, 

Then we have to prove that the charges Xj are positive measures. From the definition 
(n:^ of $ we have 

/i = (log |<l>| + 2arg$)', hQ-hj=(\og ^ +iarg^^) . 
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Recalling (12.401) and fl2.42p we have (see the notation in fl2.45p ) 



Aj = Tqj ^ |$o| = on Aj = 7j. 

This relation clearly holds locally (in a neighborhood of the endpoints of A^) and since 
Definition 12.61 and 12.71 of the case I and since the arcs jj do not intersect the cuts of Dl, it 
holds also globally on the whole of 7j. Therefore 

^;(« = 2l^"^l||- ?^^(=r„.=7,). (3.3) 
We now recall the parametrization of F = Fq,! U Fo,2 U Fi 2 see f l2.32p 

T:={z: J{iy,z) = 0, z/ G [-2,2]}. 

Since the polynomial J(i/, z) has the form (13. 2p . the parameter u is related to the argument 
of $o/$j by 



cos(arg— ^ 2' 

Hence, when ^ goes along •ja^ from aj to Pj = 7^^. (— 2) = 7f,^. (— 2), i.e., when u changes 
from 2 to —2, the argument of $o/$j monotonically changes from to tt and when ^ 
continues to go from j3j to bj, the argument of ^o/^j monotonically changes from tt to 
271 (a change of the argument here from vr to would contradict the total variation of 
the argument, which follows from the argument principle applied to the function $o/'^j)- 
Thus, the argument in (13. 3p is a strictly monotone function and therefore, with the proper 
orientation of 7j, 

a;.(O>0, eeA„ J = 1,2. 

The behavior of A^- at the endpoints of A follows from the fact that the function h has 
poles (see (11.220 ) at the points ai, bi, 02, b2. This proves Theorem 12. 10[ 

The first part of Proposition 12.91 follows from the Cauchy-Riemann relations and (13.30 



^ 1 

on 



$0 



The fact that for some A G I the set fio,i is nonemty can be verified by checking some 
particular examples of A, see Figure 12.61 second and third picture. □ 

3.1.3 Proof of the geometric results for the case II 

The following statements (Proposition 12.131 and Theorem 12.141 are related to the geomet- 
rical case II. 

Proof of Proposition \2.13[ For the geometrical case II (unlike for the case I) the analytic 
arcs 7i and 72 (i.e., the trajectories forming the contour F in (I2.32p ) do not coincide 
with the piecewise analytic arcs Ai and A2 (i.e., the cuts which make the /i and /2 
holomorphic) and they do not coincide with the cuts 5o,i ^-nd 5o,2 in (11.270 of the sheets 
of 9^ 

7j7^Aj> 7i7^Vj, J = 1,2. 
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Nevertheless, locally (in a neighborhood of the branch points of 9^) and therefore globally 
(as long as the trajectories going from the branch points do not touch other cuts of DX) 
we have, with the notation (I2.48p and fl2.45p . 

Aj := -fj n Aj C Toj ^ |$o| = on Aj, j = 1, 2. (3.4) 

Then the trajectory 71 passes through the cut 60^2 and the branch $0 goes over to $2- 
Thus we have 

71 n i?2 = -£^2 c ri_2. 

Analogously, 

72 n ^1 = ^1 c ri,2, 

i.e., 

|<l>i| = |$2| on El U E2. (3.5) 

It remains to determine which branches of $ have the same modulus on Ai^2- The curve 
7q, encloses one pair of the branch points and with the convention (12.561) this is {02, 62}- 
Indeed, it cannot enclose one or three branch points because that would contradict the 
compactness of 9^, and it cannot enclose four branch points because that would contra- 
dict the maximum principle. The fact that it encloses the pair {aj,bj} follows from the 
definition of the case I. Thus, we can join infinity and a point of Ai by a path without 
crossing 7q, \ Ai 2, and at the same time any path joining infinity and a point of A2 has 
to cross 7q, \ Ai 2. Taking into account that in a neighborhood of infinity the branch $0 
has the smallest modulus and (13. 4p . we conclude that 

7a \ ^1,2 C ri,2, 

and moreover we have that the neighborhood of infinity bounded by the parts of F belongs 
to no,i,2- Also, the connected domain bounded by the parts of F and containing {02,^2} 
belongs to ^0,2,1- Then the trajectory 7^ passes through the cut 5o,2 and the branch $2 
goes over to $o- Thus we have 

7« C Ai,2 = Ai,2 C Fi,o. (3.6) 

To complete the proof of the proposition we notice that (13.50 and (13. 6p imply (see the 
notation ([239])) 

G C ^^0,1,2- 

□ 



Proof of Theorem \2.14 



1. From Proposition 12.131 it follows that 



^0,1,2 U ^^0,2,: 



Hence the modulus |$o| is the minimal among all the moduli of all solutions of the 
equation for $. For a solution of the equation for $ one has 

log |$| G Harm(C \ {oi, 61, 02, 62}) 
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hence the function log I'I'ol is a superharmonic function in C \ {oi, 61, ^2}; and by 
the Riesz decomposition theorem, see e.g. [631 Theorem II.3.1], 

log|$o(^)| = K"(z)-log|co|, (3.7) 

see f l2.3ip . where A is a positive measure. The total mass of this measure is equal to 
2 due to fl2.3ip and this measure is supported on Aq = Ai U A2 because log |$o| ^ 
Harm(C\ A). The relation fl3.7p can also be obtained by taking the real part of the 
primitives of both sides of the Cauchy integral formula 



%{z) 1 f 



$0(2;) 27rz ^ 



Since ^q/^o = h and taking into account the pole type singularities of h at the 
branch points, we obtain part 1) of the theorem. 

2. Since the jump of the function Hq over Ai U A2 produces the positive measure A, the 
jump of the function hi over Ai and the jump of the function /i2 over A2 produce 
measures with negative sign, i.e., — Ai and — A2. Repeating the arguments of the 
proof of Theorem 12.101 we see that the jump function h2 over A2 U Ei produces a 
measure with negative sign 

-/ii - A2, supp(/ii) = El, supp(A2) = A2, 

with total mass equal to —1 

-IahI - IA2I = -1- 

The measure /ii with positive sign can also be obtained by means of the jump of 
the function hi over Ei. This proves part 2) of the theorem. 

3. Since the total mass of the charge produced by the jump function hi over Ai U Ei 
is equal to —1, see (11.211) 

lAiil - |Ai| = -1, 
we obtain the last part of the theorem 

□ 



3.1.4 Proof of the geometric results for the case III 

The following statements (Proposition I2.15[ 12.171 and Theorem 12.180 are related to the 
geometrical case III. Proposition 12.151 is just a more explicit form (relating to III) of 
the general result from Proposition 12. 4[ The equation (I2.6ip - (l2.62p in Proposition 12.151 
is well known (see [S]; the equivalent form of the algebraic equation for the Riemann 
surface of the case III goes back to Nuttall [S21 [12] )• The proofs of Theorem 12.181 and 
Proposition 12.171 are a repetition of the proofs of the corresponding results for the case I 
in Theorem 12.101 and Proposition 12.91 A minor change in Theorem 12.181 is the behavior 
of the weight A'^ at the point b*. This change is due to the regularity at b* of the function 
h given by fl2.58p . unlike in the case A G I. Another minor change in Proposition 12.171 is 
that the domain Qi^q is always present for case III; for the case I it exists just for some 
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A G I. To prove this fact we notice that all trajectories 7^. , j = 1,2,3 belong to Fq,!, 
see the notation in (12.451) . Indeed, the multiplicity of the root in (12.631) and (13.21) gives 
that all 7^*^ belong to the same r^^^, but 7^!'' C Fq,!, see Definition 12.161 -2. Thus, from 
the outside the Jordan curve J 

Jc7i'^U7lf^cro,i, b*eJ, 

is a part of the boundary of ^0,1,2 and inside it contains a tleast a part of 72 C ro,2, 
otherwise it would be in contradiction with the maximum principle, and therefore ^1,0,2 
is always inside J (see Figure [2A01) . 



3.1.5 Proof of the geometric results for the cases IV and V 

As we already mentioned. Theorem 12.211 gives a complete classification of the geometry 
for the problem with three branch points 



Proof of Theorem \2.21\ and Propositions \2.2^ \2A 



1. According to the definition (I2.68P of the domains Dj, j = 1, 2, 3, if a point a belongs 
to the boundary between the domains Di and D2, then the trajectory '■ v € 

(—2,2]} has to pass through the point 6 = (see Figure [27101) . To find this set of 
points a we take the explicit form of the equation ( 12.33p . substituting in ( 12.34^ the 
explicit expressions (12.621) : 

J(z/, 2) = 

and set z = 0. As a result we have 

(z/ - 2){a^ + 4(1 - v)a^ + (22 - %v)a^ + 4(1 - v)a + 1)^ 

.j\y^z) — 



16a2(a + 1)4 

If we set the numerator equal to zero, then we obtain the algebraic parametrization 
of Di. 

2. The boundary between D2 and 1^3 is formed by the critical triples (see Defini- 
tion [223]). To find the set of points a such that the triple { — 1,0, a} is critical, we 
can proceed in two ways. The first way is just to substitute 



9(a2 + a + 1) 

from (12.591) in the explicit form of the equation (I2.33p . As a result we get 

11664a4(a2 + a + l)6(a + l)4 ^' 

where P{a, v) is the same polynomial as in part 2) of Theorem 12.41 Thus we ob- 
tained the characterization of the critical triples for which for some v G (—2, 2] the 
trajectories of F pass through h* . The second way is to follow Definition 12.231 and 
to find out when the algebraic curve J{i>,z) has branch points for u G (—2,2). For 
this purpose we compute the discriminant of J(z/, z) and as a result we find 

Vj = C{u - 2f{u + 2)3p(a, z/)2p(a, ^/)^ 
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where C is a constant and P{a, u) is a polynomial for which the zeros do not produce 
branch points (because P appears with a square). The zeros of the polynomial 
P{a, v) are the branch points and the local analysis {zj — ~ (z/ — shows 
that the trajectories at these branch points behave like 

Z = Co + C2(Z/ - Z/q) ± C3(Z/ - Z/o)^''^ 

i.e., at these branch points the trajectories meet each other at a zero angle. This 
proves Theorem 12.211 

3. Now we complete the proof of Proposition 12.241 (part A of this proposition has 
just been proven). It remains to prove that for a strictly non-acceptable triple 
the analytic arcs 71, 72, and 73 intersect at one point c. First we notice that 
these arcs cannot join the point h* with the points 01,02,6 without intersection, 
because otherwise we get a contradiction with the irreducibility of $. Indeed, $ 
has a quadratic branch point at h* and therefore along all three trajectories 7^!'' 
{i = 1,2,3) we have the equality of the modulus of the same two branches and 
these trajectories therefore give equality of the values of these branches to all other 
branch points ai, 61, b. The same reason shows that 71, 72, and 73 cannot have only 
pairwise intersections. Thus for this triple we must have a set of points Z where the 
three arcs intersect and where all three branches have the same modulus 

Z:={z:|<l>o(^)| = |$i(^)| = |$2(^)|}. (3.8) 

If we take the normalization (12.301) into account, then we have for these points the 
prescribed values of $ 

Z:={z: {$o(^),$i(^),$2(^)} = {l,e2-^/^e^'^^/3}}. 

From the equation fl2.6ip for $ it then follows that there exist at most two such 
points. For an acceptable triple we may have all possibilities; see Figure 12.101 where 
Z = in the first picture, Z has one point in the second picture, and Z has two 
points in the third picture. For a non-acceptable triple we have that Z contains 
one point; as we already explained Z = contradicts the irreducibility of $ and 
two points in Z leads to a contradiction with the maximum principle. This proves 
Proposition I2.24[ 

4. Proposition 12.221 follows immediately from Definition 12.201 {b G f2o,i) and the fact 
that b G ro,2- 

□ 



The last set of statements (Proposition 12.241 and Theorem I2.26P is related to the geo- 
metrical case IV. Their proofs go along the same lines as the proofs of the corresponding 
statements for the geometrical case II. One minor change should be taken into account, 
namely the regular behavior of the function h at the branch point b*. It causes a change 
of the local behavior of the measures near this point. 



60 



3.2 Proof of the vector equilibrium properties for the weak Um- 
its of the Hermite-Pade approximants 



The weak asymptotic formulas and the convergence theorems for the Hermite-Pade ap- 
proximants (Theorems 12.281 12.291 and 12.301) are direct consequences of the strong asymp- 
totic formulas (Theorems 12.351 and I2.36P which will be proven in the next section. In this 
section we concentrate on the verification of the vector potential problem related with the 
weak asymptotics. 

Proof of Theorem \2.31[ In order to prove the theorem we have to verify the statements 
of the theorem for each geometrical case of the position of the branch points A := 
{ai, 6i; a2, &2}- We shall try to do this in a unique way. The proof will be given in 
the following way. 

1. First, we consider the real part of the Abelian integral (11.181) . (I1.27p 

g{i) ■= Re ^(0, e e 9^, (3.9) 

where 9^ is an arbitrary algebraic Riemann surface of order 3, and we introduce 
universal global branches for g := {(7o; fi'i? 5'2}- These global branches of g define a 
sheet structure for 9^. This universal sheet structure for 9^ goes back to Nuttall [62] 
and is different from what we use for the definition of our geometrical cases. 

2. Second, we define two universal measures A and M supported on the new cuts of 
9^. These measures have total mass |A| = 2 and |M| = 1 and they satisfy the vector 
equilibrium property (I1.12p with the matrix of interaction 




(3.10) 



The vector potential problem (11.120 with (I3.10p was introduced by Nikishin [571 [58] . 

3. Finally, we perform a reglueing of the Riemann surface and make a correspon- 
dence between the universal sheet structure of 9^ and the specific sheet structure 
{9^0, 9^1, 9^2} which we have defined for each geometrical case. In practice this pro- 
cedure gives us the relation (by means of some balayages) between the measures 
(A, M) and (A, /i) and transforms the universal vector potential problem with the 
Nikishin matrix of interaction (I3.10p to the potential problem (12.830 in the theorem. 

We emphasize that the potential problem of Theorem 12.91 also has a universal character, 
i.e., it does not depend on the geometrical case under consideration. Although the problem 
(I2.83P is more sophisticated in comparison with the potential problem (11.120 with (I3.10p . 
this problem has an advantage because it is formulated in terms of the functions /i,/2 
which we are approximating (i.e., in terms of the branch points and the cuts which make 
the functions holomorphic) and not in terms of 9^. We also mention that the procedure 
to construct the proper 9^ starting from the functions /i,/2 is the most difficult step in 
proving the Nuttall conjectures (see [621 Section 3.4]) and this was in general not solved 
before. 
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We start the proof with an arbitrary algebraic 9\ of the third order and consider 
the real part of the Abelian integral (11.271) . which is a single- valued (up to an additive 
constant) function on 9^, with the only singularities at the points at infinity 



^o(o = -2iogiei + co + ---, e^oo(o), 

9j{0 = log lel + c, + ■ ■ ■ , e ^ oo(^-), J = 1, 2, 



(3.11) 



and 

This local definition of the branches of g{^) can be continued on C as follows 

WzeC: go{z) < gi{z) < g^iz). (3.12) 

From (13.111) and the maximum principle for harmonic functions we have 

go{z) + gi{z) + g2{z) = Co + ci + C2 = 0, Vz G C, (3.13) 

and by (13.131) we have then fixed the additive constant for g. Thus we have two piecewise 
analytic curves 

Ta:= {z eC: go{z) = gi{z)} 



Tm := {z eC: gi{z) = g2{z)} 
which define the Nuttall sheet structure of IK 



(3.14) 



mo:=C\TA, ^i:=C\{T^UAm), 9^2:=C\rM. (3.15) 
We shall denote the Riemann surface with the sheet structure (I3.15P by 



2 

3=0 

Then, like in the proof of Theorem 12.141 we have by (13.121) that g^ is a superharmonic 
function in C and a harmonic function in C \ Fa, and g2 is a subharmonic function in C 
and a harmonic function in C\Tm- Therefore, by the Riesz decomposition theorem there 
exist positive measures A and M such that 



^ ^ y I —t— I r\ >Nlllllll/ll I \ 

(3.16) 



go{z) = V^{z) + Co, supp(A) = Fa 



^g2iz) = -V^\z) + C2, supp(A/r) = Tm. 
From (13.111) we have for the total mass of A and M 

|A| = 2, |M| = 1, 
and (I3.13p - (l3.14p gives the Nikishin vector equilibrium relations 

2V^{z) - V^z) = -2co - C2, on F/ 



(3 17) 

-V^{z) + 2V^'iz) = 2c2 + Co, on Tm- 
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We also note that, in the notation we used earher, we have for \z\ large enough, 

g,{z) =\og\^,{z)l (3.18) 

however, the domains where the branches are single-valued (which we defined earlier) 
are different from the domains where the branches are single-valued. We define the sec- 
ond set of global branches for the algebraic functions h = {/lo, ^i, ^2} and $ = {$0, $1, $2} 
in accordance with the Nuttall sheet structure (13.151) : 

ho, $0 G H{C \ Ta), $1 G H{C \ (Fa U Tm)), $2 e H{C \ Tm). 
Then we have 

dA{z) = ^pio+iz) -ho^{z)] dz, zeT^, 

/- ~ \ ■ (3.19) 

-dM{z) = ^i[h2+iz) - h2-{z)j dz, z G Tm. 

Now we perform the following reglueing procedure to transform the Nuttall Riemann 
surface 9^ (I3.15P into the Riemann surfaces with sheet structure defined in Subsections 
I2.2H2.41 for which we use the notation IH. We note that 

Ta U Tm = r, 

where T is defined in (I2.32p - (l2.34p and does not depend on the sheet structure of IH. The 
union of the cuts of the sheets {IHo, ^1, 9^2} coincides with F: 

2 

[J ^i,k = r, 

£,fc=0 

and at the same time the union of the cuts of {9^o, 9^2} occupies just a subset of F: 

2 

U Se,k C F. 
e,k=o 

To transform into 91 we consider the regions ^j,k,e, j ^ ^ pairwise, j, k,i = 0, 1, 2, 
see (12.441) . The contour F divides C into these regions. Then we take the components 
of flj,k,e with indices {j,k,i) 7^ (0,1,2), lift them to TT~^{Qj^k,e) on 9^, cut 91 along the 
boundaries d{'n-^^{Qj^k,e)), and interchange these pieces of 91 by reglueing them such that 

T^J^i^j,k,i) ^ T^o^i^j,k,i) 
^k^i^j,k,e) ^ T^i^i^j,k,e) 
7rg^{Qj^k,i) ^ 'n'2^{^j,k,i)- 

As a result we obtain 9^ with Nuttall's structure of the sheets. During this procedure some 
extra cuts have appeared. The transformation of 9^ to 91 works in the reverse order and 
some cuts may disappear. Since 00 G ^^0,1,2 for all 9^ under consideration, this procedure 
provides balayages of some parts of the measures A and /i to parts of the measures A and 
M. This allows to transform the potential problem (I3.17P to the potential problem (I2.83p . 

We now apply the transformation 9^ ^ 9^ to prove the theorem. We start with the 
cases I and III. Of course, for these cases it is easy to prove the theorem directly (see 
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ro,i. 

. Oi 


^0,2,1 


ai hi 

^0,1,2 


{ t y "^X / 

O ^1.2 
"1,0,2 



Figure 3.1: The partition of C by F (a generic case for I and III) 



Remark \2.32\ parts 1 and 3). However, for methodological reasons we use here the more 
sophisticated general procedure. We consider a generic case of the classes I and III (see 
Figure [2T6l second picture, and Figure [2rT0| second picture). For this case all the regions 
^0,1,2, ^0,2,1, and f2i,o,2 are present. It follows from Propositions 12.91 and 12.171 that the 
contours Fa and Tm are (see Figure [XT]) 



ro,i U 7c,f,2, Tm := Fi,2 U 7, 



c,a2- 



Thus, the sheets {9^o; 5^i; ^^2} are as shown in Figure [3^ If we re-glue the sheets IHo 
and yii along 7r~^(Fo,i \ 701,61), interchange the domains 7!-Q^{Qifi2) and 7rj"^(f2i^o,2), and 
then glue D^o and IKi along 7r^^(Fo,i \7ai,6i), and and IH2 along 7r~^(7a2^c), we see that 
^~^(7a2,c) disappears from the first sheet and moves to the zero sheet. If we do the same 
with the sheets IHi and IH2 and the domains vrjf ^(^0,2,1) and 7r^^(f2o,2,i)5 then we obtain 
the Riemann surface 9^ from Figure 12.51 

Decomposing the measure A we have, due to fl3.19p and the structure of IH, 

A = Xi + XcM + f^o,i, (3.20) 

where 

Xc,b2 •= X 



ro,i\Ai 



, fJ'0,1 ■= A 

7c,62 

and Ai, A2 are from Theorem 12.101 Analogously, we have 

M = A,,a2 + /ii,2, (3.21) 

where 



Ac,a2 •— -^2 



/ii,2 := M 



1,2 



The equality 

go = log|$o|, 

see ( I3.18p . is valid in the neighborhood of infinity bounded by F. We have 

log|<l>o| = y^-log|Co| 

due to the definition of A in fl2.46p and the relations between h and $, see (12.270 . (I2.30p . 
and (I2.3ip . Therefore in the neighbhorhood of infinity bounded by F, we have 

yAi ^ -^A,,,, ^ y^o,i + = yXi + _ log I Co I, 
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Figure 3.2: The Nuttall sheet structure of (the case I) 



and we conclude that 

Co = -log|Co|, 

and the measure /io,i is the balayage of the measure Ac,a2 from 7c,a2 to Fo^i \ Ai, i.e., 

yi^oA = YX.,a^^ Yo^i \ Ai and outside. (3.22) 

Analogously we have in the neighborhood of infinity bounded by Tm 

§2 = log|$2| 

and 

log|$2| = -K^^ - logical. 
Therefore in this domain we have 

yXca, ^ y 1^1,2 ^ yX2^ C2 = -loglCsl, 

and we obtain another balayage relation 

^ v^^-''^, on ri,2 and outside. (3.23) 

In addition to the balayage relations fl3.22l) - fl3.23p . which are valid outside Fq,! \ Ai and 
Fi^2, we need relations between the potentials inside these curves. Let us consider the first 
relation of the Nikishin equilibrium problem fl3.17p on the curve Fo,i \ Ai. If we substitute 
the decompositions fl3.20p - fl3.2ip into fl3.17p and use the balayages fl3.22l) - fl3.23l) we obtain 

2V^^ + l/^o-i + = -2co - C2, on Fo,i \ Ai. (3.24) 
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The left hand side of (13.241) is a harmonic function inside the domain bounded by ro,i\ Ai. 
Hence by the maximum principle for harmonic functions, the relation fl3.24p is valid in 
this domain. In the same way we obtain from the second relation in (13.171) that 

+ yK^,c _ y\i ^ + Co, on ri,2 and inside. (3.25) 

Now it is easy, for the cases under consideration, to derive the equilibrium relations of 
Theorem 12.311 from the equilibrium relations ( 13.17^ using ( ]3.22p -( !3.25p . The first relation 
of (IHTTl) considered on Ai gives us, using ([320]) -([321]) and (Km - (K23\f . 

2V^' + V^^ = -2co- C2=: Ki, on Ai. 

The first relation of ( 13.171) on 7^63 (see Figure [3TTi) gives us, using also ( 13.251) . 

2yA2 ^ yx^ ^ ^2c2 + Co) = -3(co + C2) =: ^2, on 7^,6^. 

We get the same relation on ^a2,c from the second relation of ( 13.171) . from ( 13.231) and with 
the help of ( I3.24p . which is also valid on 7^2,0 as we have seen. Thus, from the general 
Nikishin potential problem (I3.17P we obtain the Angelesco potential problem (11.120 . (11.150 

2^Ai + 1/^2 = Ki, on Ai, 
yAi ^ 2\/^2 = K2, on A2. 

Therefore for the cases I, III, and V (see Remark 12. 32[ 1-4) the equilibrium relations are 
verified. We emphasize the importance of the balayages ( ]3.22p -( !3.23p for the description 
of the weak limits of the extra interpolation points in the case V (see Remark 12.321 -4). 

Now we apply the same approach to verify the equilibrium relations ( 12.70 for the re- 
maining geometrical cases II and IV. For the case II Proposition 12.131 gives (see Figure [2^8] 
and ([238]), ^M) 

Fa = A1UA2, 

Tm = S1UE2UE0, ^o:=7a\Ai,2, 
and the corresponding measures (13.161) and (13.191) are given by 
A := A1 + A2, 



M := /ii + /i2 + Aio, fi2 ■= M 
where the measures 



, /io := M 

E2 



Eo 



Ai = A 



A2 = A 



A2 



are defined in Theorem 12.141 The reglueing of the Riemann surface !K gives the balayage 
relation 

^ y,,2 ^ yi^o^ on E2 U Eo and outside. (3.26) 

The second equilibrium relation of (I3.17p . the balayage (13.260 and the maximum principle 
give 

2V''' + + - V^' = 2c2 + Co, on E2 U Eq and inside. (3.27) 
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The first equilibrium relation in fl3.17p on Ai and (13.261) give the first equilibrium relation 
in fl2:83|l 

+ 2V'^^ - V^"' - V^^ - V^'' = 2V^' + V'^' - V^^ = -2co - C2 := ki, on Ai. 



If we take the first relation in fl3.17p on A2 and use (13.271) . then we find the second relation 
in fl2:83|l 



2^' + 2V'^^ - V^' - V^^ - V^'° = V^' + 2V'^' + =Ki + 2c2 + Cq := ^2, 



on A5 



Finally, the second relation in (13 . 1 71) on E = EiU E2 and the balayage (I3.26P lead to the 
third relation of fl^TS^ 



2\/^i + 2\/^2 + 2l^«' - - V^'' = 2^' + V^^ - V^' = 2c2 + cq = K2 - ^i, 



on E. 



Thus the equilibrium relations (12.831) for the case II are verified. 

If we verify the equilibrium relations for the last case IV, then we will also get the dual 
equilibrium problem (see II and III of Theorem 12. 31 p . For the case IV Proposition 12.251 
(see dam]) and (l2T5|) ) gives 



AiUA2 = AiUAi,2UA2, 
£-0 U El U E2- 



The Nuttall sheet structure of 9^ is shown in Figure 13.31 It defines the global branches 
(I3.16P for the real part of the Abelian integral (13.110 



go = v 



+ Co, 



= _\//»l+/»2+/i0 

where the measures 

Ai = A _ , A, = A 



Co - C2, 



~ 5 

A2 



Ai,2 = A 



^l2 = M 



E2 



are the same as in Theorem 12.261 and 



/ii := M 



El 



/xo :=M 



Equating 



qq = gi, on Fa, 
gi = g2, on Tm, 



gives the Nikishin equilibrium (13.71) . Observe that (13.71) has just two equilibrium relations 
since the 0-sheet and the second sheet do not intersect. 
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Figure 3.3: The Nuttall sheet structure of 9^ (the case IV) 

In order to obtain (12.831) we re-glue the Riemann surface shown in Figure 13.31 We 
interchange the domains bounded by Eq U E2 on the sheets 1 and 2 and as a result the 
cut along Eq U E2 between the sheets 1 and 2 disappears and we arrive at the Riemann 
surface 9^ shown in Figure 13.41 

The new sheet structure defines new global branches for the real part of the Abelian 
integral (13.111) 

f 90 = + Co, 

g=: h, = + yMi _ Co - C2, 

[^2 = -F^^ -F^^ +C2. 

Equating ^0 = 9i on Ai, where Ai is the cut joining the sheet and 1, we have 

^ y\2 ^ yXi,2 ^ Co = -V~^^ - V^'-^ + V^"' - Co - C2, on Ai, 
which gives the first relation in (12.831) : 

2V^i + _ = _2co - C2 := Ki, onAi. 
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Figure 3.4: Reglueing the Nuttall Riemann surface (the case IV) along Eq U E2 

Equating ^0 = 92 on A2 gives the second relation in fl2.83p 

yx^ ^ 2V'^, ^ y^., =c2_co =, on A2. 

Finally, equating gi = g2 on Ei gives the third relation in fl2.83p : 

_YXi ^ yX2 ^ 2V^i = 2C2 + Co = - ^i, OU ^i. 

Thus the equilibrium relations (12.831) are verified for the case IV. 

To derive the dual equilibrium problem we re-glue the Riemann surface shown in 
Figure [331 by interchanging the domains bounded by £'0 U Ei on the sheets 1 and 2, and 
as a result the cut along Eq U Ei between the sheets 1 and 2 disappears and we arrive at 
the Riemann surface shown in Figure I2.14[ The new sheet structure defines new global 
branches for the real part of the Abelian integral (13. lip : 

[ 90 = + V^' + Co, 
9=- lgi = -v'^' - V^^ - Co - C2, 
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Equating = on Ai (see Figure gives the first relation for tlie dual equilibrium 
problem 

2V'^' +V^' + V''^ = Ki, onAi. 
Equating = §2 on A2 gives the second relation 

yXi ^ 2V^^ - V^^ = K2, on A2, 

and Qi = §2 on E2 gives 

yXl _ y\2 ^ 2V^^ =Ki- K2, OU E2. 

Thus the equilibrium problem for the case IV (see II and III of Theorem I2.3ip is verified. 
The dual equilibrium problem can be verified in a similar way for the other geometrical 
cases. 

To conclude the proof of Theorem 12.311 we verify the symmetry property (12. 84^ . We 
show how to do this for the geometrical case I; one can use a similar reasoning for the 
other cases. The equilibrium relations for the case I have the form 

Ui ■= 2V^' + V^^ = Ki, onAi, 
U2 := V^' + 2V^' = K2, on A2. 



For the harmonic conjugate Uj one has 



dn dr 

hence one has to verify for Uj = Uj + iUj that 



0, on Aj, j = l,2, 



-^j (^) = T^^i (^) ' ^ e A, \ {a„ } , J = 1,2. 



This relation holds because the function Uj is analytic on A^ \ {aj,bj} (because Uj = 
log$o/'^'i + const), and therefore its derivative does not depend on the direction. The 
theorem is now proved. □ 



4 Asymptotic analysis of the matrix Riemann-Hilbert 
problem for the Hermite-Pade polynomials 

In this section we prove the results about the strong asymptotics of the Hermite-Pade 
polynomials, i.e.. Theorem 12.351 and Theorem 12.361 We demonstrate here the detailed 
proof for the case II. The proof for the case I is just a simplified version of the proof for 
the case II and the most delightful pieces of the proof degenerate. The same is true for the 
cases III and V which are both a simplified version of the proof for the case IV, and the 
proof for the case IV is just a repetition of the proof for the case II with minor differences 
in the local analysis of the asymptotics in the neighborhood of the branch point b*. For 
these reasons we present the proof for the generic case II in detail and give a sketch of 
the proofs for the other cases. 
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For the proof we use the steepest descent method of Deift and Zhou (see [21] )• The 
method consists of several consecutive transformations 

of the matrix Riemann-Hilbert problem (I2.13p for the matrix function Y to the Riemann- 
Hilbert problem 

' Z e /73x3(c\S), 

3Z± G C{t) ■.Z+ = Z^[I + 0(l/n)) , uniformly on S as n ^ oo, (4.1) 

{Z{z) = I + @{l/z), z^oo. 

It is known that for large n, the problem (14. ip has a solution which satisfies 

Z = I + 0(l/n), uniformly in C \ S as n — oo. 

Thus, after making the inverse transformations from Z to Y we obtain the existence of 
the solution Y of the Riemann-Hilbert problem (12.13P for large n and asymptotics for 
n — s> oo for its components. 

4.1 Normalization of the Riemann-Hilbert problem at infinity 
and decomposition of the jumps 

The goal of the first step is to transform the Riemann-Hilbert problem ( 12.13p -( l2.17p such 
that the solution of the new problem has the same normalization at infinity as the solution 
of ( 1411) . For this purpose we use the algebraic function $ = {$o, $1, $2} in (l2.27p -( l2y3Tl) . 
We denote 

/<1>[J \ /Q \ 

^ := <1>^ , C := , (4.2) 

\0 0<1>^/ \0 C^J 

where Cq, Ci, C2 are the inverses of the leading coefficients of the expansion of $0) "^i; '^'2 
near infinity, see (I2.3ip . If we set 

Z := CYS, (4.3) 

then, because of d^TTS]) and (lOTD 

Z(;2)=/ + 0^, ^-^00, (4.4) 

and the piecewise analytic matrix function Z satisfies the following jump condition: 

Z+ = Z_J, on Ai U Ai,2 U A2 U ^1 =: S, (4.5) 

(see Figures 12.81 and 14. ip . and for the jump matrix J we have from (12.140 that J : = 
SZ^WS+ so that 

'J(wi,W2), on Ai,2, 
J{wi, 0), on Ai, 



J :-- 



(4.6) 

J{0,W2), on A2, 
tj(0,0), on El, 
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where we set 



J{Wi,W2) 



v 



0- 






Wi 



0- 

$^ 





2+ \ 

0- 





(4.7) 



Note that, in comparison with Y, the RH problem for Z has a new contour E'l where Z 
is discontinuous, due to the jumps in the functions $j. 




Figure 4.1: Jump contours for the Riemann-Hilbert problem for Z and Z (normalization 
and the global lens) 

Now we analyze the jump matrices (I4.6l) - fl4.7p . We start with the usual decomposition 
of the jumps on Ai and A2, which are commonly used in the analysis of 2 x 2 matrix- valued 
Riemann-Hilbert problems. If we define 



Di-.-- 



D 



2 •- 



/ 1 





o\ 


n 


1 





V 





V 


/ 1 





o\ 





1 





V$>2 








Wo 



/ 





Wi 






1 
















Wi 




1/ 


\ 


















( 





1 






(4i 



1 

W2 



— 00 
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then we have, since $o± = *^i=f Ai, 



J = J{wi,0) 



(in 

^1 



\ 








Wi 

n 







1/ 



on Ai, 



(4.9) 



and similarly, since $o± = *^'2=f ^2, 

J = J(0, W2) = D2+W2D2- 



on A 



2- 



(4.10) 



The following decompositions of the jumps on Ai 2 and Ei were not used before in the 
analysis of Riemann-Hilbert problems. If we define 



1^2,1 : = 



\ 

$21^2 



/I 
1 
\0 1 



then we observe that 

J = J{wi,W2) = D2,i-J{wi,0)D2}+, on Ai,2, 



(4.11) 



i.e., after the decomposition (14.111) the jump has a block structure and J{wi, 0) can then 
be decomposed again in the usual way as in (14. 9p . since $o± = '^*i=f holds also on Ai 2. If 
we define 



(\ 





V 




1 



0\ 



1 



/I 




I 

then we have on E\ (since $i± = $2^= 

^10 \ 




— 

W2 



\ 

W2 





(4.12) 



J= J(0,0) 



Since 



1 



v 



$2^2 







/l \ 



.0 1 








-- 



1 
\0 



_$>2_ 

1 



(\ 








1 







^2:1+- 



V 2<l>>2 / 




on 



we have 



J= J(0,0) =Di,2-W^l,2 



/I 









1 







V 2$>2 
73 



1 



-^^2,1+5 



on i?!. 



and we obtain 

J = J(0, 0) = Di,2-W^i,2/^i,2+/^2:i+, on E,. (4.13) 
We also point out the commutation relations 

D1D2 = D,D,,2 = Di,2Di. (4.14) 

4.2 Opening a global lens (in preparation of opening local lenses) 

The goal of the second step is to transform the Riemann-Hilbert problem f l4.4l) -( HI71) for 
the function Z in (14. 3 P to a Riemann-Hilbert problem with new jumps which can be 
decomposed in the usual way by (I4.9I) - (I4.10I) to jumps which do not depend on n or to 
jumps which tend to the identity matrix as n ^ 00. 

We fix a neighborhood Ocj of the point Cj (j = 1, 2) (see Figure \AA\} and denote 

A^-) := AanO,^.. 

We join the endpoints (different from Cj) of Ag'^^'', which we denote by c^, by an arc E2 
such that 

C n2,i, (4.15) 
where we used the notation (12.45^ . We denote by the domain bounded by 

dG+ := El U Af^^ U A^'^^ U E^. 

(see Figures 12.91 and 14. ip . We define 

Z := < _ ' 4.16 

\Z, inC\G+. ^ ^ 

This piecewise analytic function Z has a jump 

Z+ = Z^J, on S U E^ =: S+ (4.17) 
where for J on S+ \ (Ai,2 U U E+ U A^'^^ U A^"'^), see M . 

J = J={T."'l onA,U{A2\[jm- (4.18) 

On other parts of S the jump is changed to 

j^(j{0,W2)D2,u onA^^UA^), 

[D^l on E+. ^ ■ 

On El we obtain, using (14.130 . 

Z+ = Z+D2,l+ = Z_iDi^2-Wi,2Di,2+D2}+)D2,l+ = Z^Di^2-Wi,2Di,2+, 

and therefore we have 

J = Di_2_Ty, 2A,2+, on El. (4.20) 
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On Ai_2 we have, using fl4.1ip . 

Z+ = Z+D2,i+ = Z.{D2,i-J{wi,0)D2l^)D2,i+ = Z_J{wi,0), 
and therefore we have 

J = J{wi,0), onAi,2. (4.21) 

Thus the desired form of the jumps (I4.18p -( l4.2ip is indeed obtained after the second 
transformation. 

4.3 Opening local lenses 

The goal of the third step is to transform the Riemann-Hilbert problem fl4. 171) - (14.21 1) for 
the function Z from fl4.16p to a Riemann-Hilbert problem with jumps which do not depend 
on n or which tend to the identity matrix as n — > oo. We start with some notation. In 
addition to the points 

c+ := aO,^. n A2, J = 1,2, 
(see Figures 14.21 14. 3p we denote 

cj := dOc^ n El, 5+ := dO^^ H Ai,2, dj := 80,^ H Ai, j = 1, 2. 

We define a local perturbation of the arcs A2 joining the points 02 with and 62 with 
by Jordan arcs 7^^ -± and 7^^ -± 

Analogously we define (see Figure 14. 2p 

Around the piecewise analytic arc Ai U Ai^2 we now define the lens shaped domains 
and bounded by 

ar±:=AiU(A^UAt2), 

and 

Ti := T^UT{- U Ai. 
Analogously, around the analytic arc A2 we define domains bounded by 

dT^ := A2 U A^, 

and 

T2 := T+ U T^- U A2. 
Also, around Ei we define the domains T^^ 

OTi^ ■=EiUEf, 

and 
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ai 



II 




Figure 4.2: Jump contours of Z (local lenses) 

Now we can transform the Riemann-Hilbert problem fl4.17l) - fl4.21l) . We define (see Fig- 
ure SJ}, 





in T+\T2, 


ZD^, 


in T^\Te„ 


ZD^\ 


in T+\T,, 


ZD2. 


in T2" \Ti, 




in T+\Tf, 


< ZDi,2, 


in Tj;;\Tf, 


ZD^'D2, 


inT+nT^-, 


ZD^'D^\ 


inT+nT/, 


ZD,Di,2, 


inTf nT^^, 


ZD.D^^l 


inTf nT+, 




in C \ (Ti U U T^, 



The piecewise analytic matrix function Z has a jump (see fl4.17p ) 

Z+ = Z_J, on S+ U STi U dT2 U 9Tb, =: S. (4.23) 

We now describe the explicit form of the matrix function J on the different parts of S. 
From (14.201) we have on Ei \ Ti 

Z+ = Z+D^l_^ = Z_D^^2-Wi^2Di,2+D^i+ = Z^Wi^2, on \ T^. 
The same relation holds on EiHTi since Di+ = Di_ on Ei. From (HTTSD . flOT]) and (OD 
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Figure 4.3: The function Z (close-up around C\) 



we have on Ai \ (T2 U T^J 

Z+ = ZDll = Z^Di^WiDi+D^l = Z^Wi, on Ai \ (T2 U T^J. 

The remarkable fact is that this relation in fact holds on Ai fl (T2 U T^J as well. Indeed, 
on Ai n (T2 U ) we have (see Figure 14. 3p 

Z+ = Z+D^lD2+ = Z4Di^Di,2DiiDil){D,^WiD^l)D^lD2+ = Z^D^^lWiD2+. 
Note that 

D^iWiD2+ = W,, on Ai, 

and therefore 

Z+ = Z.Wu on Ai n {T2 U Te,). 

Analogously, since 

Di^2WiD^l = Wi, on Ai,2, 

we have 

Z+ = Z+D^lD^l = Z^Di^WiDi+D^lD^l 

= Z_D^_D^l{D^^2WiD^l) = Z_Wi, on Ai,2 n (T2 U T^J. 
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-1 n-1 

2 • 



In the same way we have on A2 \ Ti, see (14.101) : 

Z+ = Z+D:,l = Z_D2-W2D2+D^1 = Z^W2, 
and on A2 n Ti, see gUD and KW . 

Z+ = Z+D^'D^^ = Z_{D2W2D2+D2A+)D^lD 
Using the commutation relations (I4.14p and the observation that 

DiW2D2+D2,i+D:,lD:,l = W2, on A2, 

we arrive at 

Z+ = Z_Dl^D2-W2 = Z_W2, on A2 n Ti. 
Summarizing, we have for the jump matrix J in fl4.23p 



on A2 \ri, 



J 





on 


Af, 


D2, 


on 


A^, 




on 




< ^2,1, 


on 






on 


Ai, 


W2, 


on 


A2, 


W,,2, 


on 


El. 



(4.24) 



We also note that when we go around the point cj (j = 1,2), then the total jump of Z is 

W2W^^W{-^^Wi = I, 

i.e., the intersection points Ci and C2 are not branch points of Z. Note also that the 
matrices Wi,W2,Wi^2 do not depend on n and due to Proposition 12.131 (see Figure [231) 
and (14.151) the matrices Di, D2, D12, D2^i tend to the identity matrix I on compact sets 
of the domains of their definition in (I4.24p . Therefore we achieved the goal of the third 
t r ansf ormat ion . 



4.4 Parametrix away from the branch points 

The goal of the fourth step is to construct a solution for the model Riemann-Hilbert 
problem with jumps that do not depend on n, so that later on we can eliminate these 
jumps in (14.241) . We are looking for a function X for which 



X e iJ3x3(c\s), 

X+ = X^W, on S, 
X{z) = I + 0{l/z), z 

where S := Ai U A2 U Ei (see Figure Hi]) , and 



W :-- 



(4.25) 



00, 





on Ai, 




W2, 


on A2, 


(4.26) 


W^l,2, 


on El, 
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see the notation in fl4.8l) and fl4.12p . The function X should also satisfy the endpoint 

condition fl2:T5D - (12X71) . We used the notation S:= S \ {ai, bi, as, 62} in fOSD . 

To construct the solution of fl4.25p - fl4.26p we use the Szego functions (see subsec- 
tion [22]), i.e., the solution of the following system of scalar boundary value problems, see 

(ES]), (EH, 



1) 



2) 



3) 



FoGi/(C\(AiUA2)), 
Fi Gi/(C\(EiUAi)), 
F2G//(C\(EiUA2)), 

o 

Fo± = iFi^Ui-Wi, on Ai \ {ai, 61} =:Ai, 

O 

Fo± = iF2^uj2-W2, on A2 \ {aa, 62} =:A2, 
Fi± = F2^^^^, on El, 

{normalization condition: FqFiF2 = 1 in C, 
local behavior flCTD -2b. 



(4.27) 



Here Uj is the branch of uj'j{z) = {z — aj){z — bj), see f l2.87p . The solution of this system 
of boundary value problems exists and is given in Theorem 12.341 We define 

Fiz) := diag(Fo(2:), F,{z), ^2(2;)), F^ := F(oo), 

and transform the Riemann-Hilbert problem fl4.25l) - fl4.26p to the Riemann-Hilbert prob- 
lem for the function 

X := FZ^XF. 



(4.28) 



We have 



X+ = X_H, on S, 
X{z) = 1 + 0(1/2), 



00, 



where the jump 
due to KTnf -2. is 



H :-- 



( 







i7 = FZ^WF^ 
\ 



wx 












(\ 












/ 
\ 

<JJ2- 













V 



V 













1/ 



on Ai 



on A5 



(4.29) 



on E\. 



^1^ 
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For the solution of the Riemann-Hilbert problem fl4.28p - fl4.29p we take the Riemann 
surface 9^ from Definition 12.111 and define on 9^ the rational function Xj with one pole 
(and one zero) 



(0) 3^(1) 3^(2) 



m) = 0(1/0, e 



oo 



oo 



(0) 

U) 



Then the function 



X 



Xofi ^0,1 3;o,2 

5^1,0 5^1,1 3:i^2 

.X2,0 5:2,1 ^2,2, 



( 1 


!_ 


]_ 












[rX^^^ 


1 

A-o 


1 0.(2) 
A-o 



1,2. 
(4.30) 

(4.31) 



is the solution of the Riemann-Hilbert problem fl4.28p - fl4.29p . Indeed, the normalization 
at infinity in fl4.28p clearly holds because of fl4.30p . To verify the jump condition in 
fOHjl - fll^ we check the relation 

X+ = X_H 

O 

on the different parts of S- Substituting here (14.291) we need to check that 

/- 



' ^0,0+ 3;o,l+ Xq^2+ 
^1,0+ ^1,1+ ^1,2+ 
.X2,0+ 2:2,1+ 3:2, 2H 



\-iuJi-X2,l- 
/—iuJ2-Xo,2- 
iuJ2-Xi^2- 
\-iuJ2-X2,2- 



5o,o- 



1 

Xo.l- 



X0,2-\ 
^1,2- 
^2,2-/ 



on Ai 



3^1,1- 
^2,1- 



T-^Xi 0- 



on A 



2, 



(4.32) 



/~ UJ2— ~ ^1— ~ \ 

~ UJ2— ~ UJl— ~ 



. _ ij2^ ~ 

\X2fi- ;tt^2,2- 



UJ2 

'—^2,1-/ 



on El. 



If we substitute the expression fl4.3ip for Xk/ {k^i = 0, 1,2) into (14.320 . then we indeed 
see that fl4.32p holds identically. Thus the jump condition in fl4.28p - fl4.29p holds and X 
in fl4.3ip is the solution of the Riemann-Hilbert problem fl4.28p - fl4.29p . Finally we have 
from fl^:^ and fOO]) that 



/ Fo(oo) 



X :-- 



FooXF 



iXl 



^0, , 
(0)^1 (oo) 



iX. 



^0 ^ 
(0)F2(00) 



1 Fo(oo) 
iuji Fi 

1 _y(l) -^lM 

1 ^(1) ^2(00) 
— ^2 — F — 

LOi Fi 



1 i^o(oo) \ 



11^02 F2 
±^{2)Fi{^ 

UJ2 F2 

1 ^(2) i^2(00) 

^2 ^ 

UJ2 ^2 



(4.33) 
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is the desired solution of the model Riemann-Hilbert problem fl4.25p - fl4.26p . 



Remark 4.1. We recall that in Theorem 12.341 we presented the solution of the scalar 
boundary value problem fr2.9ip on the Riemann surface by means of the Cauchy integral 
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with a meromorphic differential on There is another way to solve the scalar boundary 
value problem fl2.9ip : see [7J for a similar approach. The idea is the following. Since 9^ 
has genus zero, it is conformally equivalent to C Let C be a conformal mapping 

cm = c, c(oo(°)) = oo. 

We consider 91 from Definition 12.111 (see Figure [2?7j) . The image of Dl is given in Figure HiH 
Analyzing the jumps around the points C(cj), see the notation fl2.88p . we see that the 
special form of the jumps (I2.89P provides the existence of a continuous solution in the 
neighborhood of the points Cj {j = 1, 2). 




Figure 4.4: Conformal mapping of to C (case II) 



4.5 Local parametrices 

The function X is not a good approximation near the branch points A := {ai, bi, a2, &2}- 
We need a local analysis near each of the branch points. The goal of the fifth step is 
to find a matrix function f/ in a neighborhood Og of each branch point e G A which 
asymptotically (as n oo) matches X on the boundary dOe and which has the same 
jumps as Z. This function is the solution of the following Riemann-Hilbert problem: if 
Cj denotes aj or bj {j = 1, 2), then we want to find a function such that 

'l) U,^eH'-%0,^\{A^UAluAj)), 

< 2) Ue^+ = Ue^J on % := O^^ H (A,- U A+ U At) : (4.34) 
3) Ue^ = (/ + &{\/n))X uniformly on dOe^ as n oo, 
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and the elements of the matrix f/g^ (z) have the same hmiting behavior as the elements of 
Z{z) as z — > Cj. 

We recall the explicit expressions of the jump matrix J, see (14 .24^ and (14.81) . 



J 





( 1 





o\ 




n 


1 


u 








\ 





1/ 




1 


Wi 


o\ 




1 










Wi 






< 


\ 





V 




/ 1 





o\ 







1 







n 













( ° 












1 







1 

V W2 





V 



on a!^\ 



on Ai, 



on A 



(±) 

2 ; 



on A2. 



We see that the jump matrices related to a branch point from A have a 2 x 2 block 
structure. The solution of the Riemann-Hilbert problem (I4.34p in the 2x2 case, with the 
local behavior of the equilibrium measure as in (I2.54p . was obtained in [17]. Thus, using 
the results in [17], it is easy to write down the explicit form of the solution for the local 
Riemann-Hilbert problem (14.341) . This solution is 



1,2, 



(4.35) 



where the matrices Ap . are 



A 



ei 



A 



62 




and, see (12. Sp . 



1/2 
1/2 



The branch of the square roots above is chosen such that 



W 



31 



on 



Aj n Oe^., G {aj,6j}, j = 1,2. 
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The matrices Eg. in (I4.35P are 

'1 i 



Ea, := ^Xdiag{wa,,Wa,^,l) y 1 j diag |^0F7I^, -^^=z, 1^ , 

1 ~' ^\ f 1 \ 
'■= 7;^dmg(wb^,w,\l) \ -I 1 diag y/nnipi, 1 

2 \0 1/ ^ v^™^ ^ 

1 ^ ' °\ ^ 1 \ 

■= 7^^diag(w)a2,l,w„2^) U 1 diag y/7m(p2, 1, , , 

2 \0 1/ ^ 

1 ~' ^\ f I \ 

Eb2 ■= o^diag(w;62,l,w;J) -2 1 diag V^V^a, 1, 

2 V 1 / ^ y^f^/ 



Here X is given by (I4.33p and we use the notation 



•= log (1^ ) • 



5 



To present the exphcit expression for the matrices Ve^ in (14.351) we follow [37] and introduce 
the matrices 

1, 

VI/,, := \nm^J',^{^) -n^,K',^{^) 


^a. := I 1 

I /,,(^) 

^fe := 1 

where la and Ka are the modified Bessel functions of order a (see 
We denote by O^e^^ the sector domain bounded by 



ao(+) := dOe, u (A,- n OeJ u (A+ n 



Analogously, oL ^ has the boundaries 



90(7) := aOe,, u (A, n OeJ u (AT n OeJ, 



and 



:= Oe, \ (0i+) u O':-^^ 



3 
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The matrices V^- in fl4.35p are 
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This gives all the ingredients for the solution fl4.35p of the local Riemann-Hilbert problem 



4.6 Final transformation. Asymptotic formulas 

We now finish the transformation of the original Riemann-Hilbert problem fl2.13l) - fl2.17l) 
to the Riemann-Hilbert problem (14.11) . The final transformation is 



ZX ^, in C \ IJj=i (^(ij U ^bj), (away from the branch points), 
ZU~^, in Oej, Cj G {aj, bj},j = 1, 2, (near the branch points). 

We have on Ai U A2 U Ei 

z+ = z^jx-^ = z^xz\x^jx-^). 

From fOSD - fOBD we have 

x^jx-^ = x^w{x.w)-\ 

so that 

Z+ = Z^, on (Ai U AaU^i) \Oei. 
Analogously, see (I4.34p . we have 

Z+ = Z_, on (Ai U A2 UEi) \0e2. 
Thus, denoting (see Figure 

S := Af U Af U A^2 U U U |J dO 



eieA 



(4.36) 
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we have 



2) Z+ = Z_J, on S, 

3) Z{z) = I + 0{l/z), z^oo. 




Figure 4.5: The contour S for the final transformation 
Since as n ^ oo 

J = / + 0(l/n), uniformly on dOe^, 



e,GA 



and for some C > 



J = I + 0(e-^"), uniformly on S \ |J dOe 



(4.37) 



we can conclude that the solution of the problem (14.371) exists for n large enough, and 

Z = / + 0(l/n), (4.38) 

uniformly for z & C 

Now we can return through the sequence of transformations (14.361) . (14.221) . (I4.16p . 

Z ^ Z ^ Z ^ Z (4.39) 

and we can conclude that for large n the solution of the Riemann-Hilbert problem (12.131) 
exists, i.e., the Hermite-Pade approximants exists and the indices are normal, and we 
get the asymptotic formulas of Theorems 12.351 and 12.361 To obtain these formulas we 
note that in the transformation (14.391) it is sufficient to keep track of the first row of the 
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matrices. We start with the asymptotics on compact sets in C \ S, where S := T,U E2, 
see (14.51) . Comparing the inverse and direct transformation in fl4.39p . we have 

Z = {I + 0{l/n))X, onC\S, 

and 

'Q$^P„ _ _ 

Z = Z= \ * * * I , in (C \ G) \ (Ai U A2) 

* * * 

'Q$^P„ CS^^rII^ C^^^R^n^ - CS^^Rll^W2/wi 
Z = ZD2,i = I * * * I 5 in G \ Ai^2- 



This gives 



Pn = 7^(1 + 0(1/^)), onC\S, (4.40) 

fli" = ^(1 + 0(1/")). onC\E, (4.41) 

-4^(l + 0(l/«)), on(C\G)\(A,UA2), 

RS> = { 7p X ^ . (4.42) 



The asymptotic formulas (I4.40l) - (l4.42p are vahd uniformly on compact subsets of the 
indicated sets. If we substitute (I4.33P in these formulas, then we arrive at the formulas of 
Theorems 12.351 and 12.361 for C \ E. 

Next we consider the asymptotics on compact subsets of S \ {ai, 61, 02, &2}- Since 
-^0,0; '^'o ^ HiC \ (Ai U X2 U A12)), the asymptotic formula (l4.4Up holds on compact 
subsets in C \ (Ai U X2 U Ai 2)) and therefore in particular on EiU E2. For the same 
reason the asymptotic formula (I4.4ip holds on compact subsets of C\E, hence in particular 
on i?2, and the second asymptotic formula in (I4.42p holds on E2. Thus the asymptotic 
behavior on E2 is verified. 

Now consider the asymptotics on Ei. On Ei we have, see (14.250 - 04.260 . 

Wi 

Xo,i+ = Xo,2-- — , $i+ = $2-, on El, 

W2 

hence we have from (I4.42p 

= + 0(1/-)) = ^^7(1 + 0(Vn)), on E,. (4.43) 



In Tbj, see (14.220 . we have 



Z 



ZD2,^D-l inTi;+\ 
ZDi,2, inTi;-\ 
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and therefore 



W2 Wi 



* 



* * 



on T, 



(+) 



El ' 



I V * 

On the other hand we have 

Z = {I + 0{l/n))X, 

hence using (14.431) we arrive at 



on T; 



(-) 



El ■ 



J 



inC\|J(a,UO,J, 



X, 



0,2+ Wi 



-^0,2- UJi ^ ^0,1- 



(l + 0(l/n)) 
(l + 0(l/n)), 



on El. 



Thus the asymptotics on Ei are also verified. 

Now consider the asymptotics on Aj (j = 1, 2). In Tj, away from ci and C2, we have 



Wi 



* 



^C^^^Pn + C^^Rl^^ Q^^^A^'^^ 



on T, 



(+) 



ZD. 



Wj 



\ 



* 
* 



on T, 



{-) 



Taking the hmiting value we obtain 



d(1) 



^(i + 0(Vn)), 



(2) 
nit 



2± 



•-^O ^2± 



(l + 0(l/n)), 



on 



A,-,j = l,2. (4.44) 



Using these asymptotics we arrive at 



^0,0+ _^ Xqj 



'-^0^0+ 



X, 



0,0- 



'-'0^0- 
Xo,o 

Q<fo/ + 



+ 



Xo,o 

'-'0^0 



(l + 0(l/r2)) 
(l + 0(l/n)) 

(l + 0(l/n)), 



on A,, J = 1,2. (4.45) 
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Finally we consider the asymptotics on Ai 2. On A^^^ we have the limiting values 



* * * 

\ * * * 



± 



From this we have that the asymptotics (14.451) for P„ and (14.441) for R^l are also valid on 

(2) 

A12. For i?n we have 

Co"$^i?S - = ^0.2(1 + 0(l/n)), in Ai,2. 

wi Cq "P^ 

Taking into account that |$o| = < \^2\ in Ai^2 (Proposition 12. 131) we obtain 

Rnl = '^^^a + mH), inA,2. (4.46) 

It remains to observe that the limiting values of the asymptotic formulas (I4.43I) - (I4.46I) 
along S at the points ci and C2 coincide. 

This proves Theorems 12.351 and 12.361 for the case II. 

Remark 4.2. We can get formulas for the local asymptotics for P„ and Rn^ {j = 1, 2) 
from (14.381) by using the transformations (I4.39P in Oej, with ej G {aj, bj}. These formulas 
are similar to the corresponding formulas from [47] and contain Bessel functions. 



4.7 Sketch of the proof for the other geometrical cases 
Case I 

For this case Ai fl A2 = 0, hence we do not need the second transformation (14.131) . i.e., 
we do not need to open a global lens. All the other steps are as in the proof for the case 
II. We just have to use for {$o; ^i; ^2} the algebraic function (I2.27p . (12.421) and we have 
to drop some parts of the proof which involve the intersection Ai fl A2. 



Case IV 

This is also a generic case. All the steps of the proof above are present for this case. The 
proof now uses the Riemann surface (I2.72p and for the normalization we use the branches 
f l2T3|) of the function 

Z := CYS, 

where we use the notation (14. 2p . The second step is to open the global lens around Ai 2 in 
the domain with boundary dG^^^ = E2U (see Figure IT6] and the notation (12.750 ) 
and we obtain 




ZDi,2, inG'+, 
Z, in C \ G+. 



Here we use the notation (14.120 . The function Z has the jump 

Z+ = Z_J, on Ai U A2 U ^2 U 
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where 

and, see fl4.8p . 



J :-- 



J := Di^2, on Ef, 
L'2-W^2^2+, on A2 := A2 U Ai,2. 



If we use the analog of fl4.13p 



then we have 



J(0, 0) = D2,i-Wi,2D2,i+D^i+, on E2, 
J := -D2,i-W^i,2-D2,i+5 on E2. 



(4.47) 



(4.48) 



ir 






















c .JE2 f^* 


ai 







Figure 4.6: The global lens (case IV) 

The third step is the opening of the local lenses around the arcs Ai, A2, E2 (see Section 
I4.3p . The decompositions fl4.47p - fl4.48p give as a result a new function Z for which the 
jumps on these arcs do not depend on n 



with 





= Z_J, 










on 


Ai 


J = W := < 


W2, 


on 


A2 




Wl,2, 


on 


E2 



and on the lenses Af , Af 2, A^, E2 



^ and E^ (excluding the points cti, 02, b, b* 
<■ 00 (due to Proposition 12. 25p . 



the jump J 

tends to the identity matrix as n 

The fourth step (the parametrix away from the branch points) is the solution of the 
Riemann-Hilbert problem for the function X with the jump W. This is just a repetition 
of Section 14.41 but exchanging Ei and E2 and with another Riemann surface 9^, namely 
fl2T2|) (see Figure EH]) . 

The fifth step (local parametrices) has a new feature which is not present in the case 
II. The solutions of the local matrix Riemann-Hilbert problems 04.340 for the function 
Ue, with e G {ai,a2,b}, are the same as for the case II (see fl4.35p ). However, around 
the branch point b* the solution is different. Because of Theorem 12.261 -2. the solution 
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is represented by means of Airy functions. The local Riemann-Hilbert problem in the 
neighborhood of b* is 

fl) f/,. Gi/(3x3)(0,*\S6*), 

2) = f/6.„ J, onS,., (4.49) 

3) Ub" = (/ + 0(l/n))X, uniformly on dOb* as n — > oo, 
where (see Figure IT71) 

Sfo. := Ob* n {E2 U U U E^), 

and, see (gSD and fICTD 

'Da.i, onE^nOfo., 

^:=<W^i,2, onEsHOfe*, (4.50) 

Di,2, onE+nOfo.. 




Figure 4.7: Local analysis around b* 

Observe the block structure of this matrix Riemann-Hilbert problem, so that we can 
use the known solution for the 2x2 case. Following [3T[ [26] we have 



where 



Ah* := diag 1, 



Ub* =Eb*Vb*Ab* 



(4.51) 



Eb* 



-v/^e^^/^X diag 1^1, (^^^ 



-1/2 / \ 1/2^ 



W2 



X I 1 1 j diag^l, (^y(/.i,2 



3n 



-1/6 



3n 



1/6^ 



90 



and 



To present we use the notation 



1 *i 

yj2 



/I \ 

Ai((f^,,)^/^) Ai (e^ (f v^,,)^/^) 

1^0 Ai'{{-^^,,f') elA^ {el{f^,,f') ) 



/I \ 



Ai((fv,,,)^^^j -e^Ai (el(f v^,,) 
^0 Ai' ((f^,,)^/^) -Ai'(ei(f^,,)^/^j ) 
where Ai is the usual Airy function and 

a ■ = diag ( 1 , e^"/^ , e"^"/^) , €3 = e^"^/^ . 

The matrix Vb* is then given by 



/l \ 

1 -1 

Vo ly 

/l o\ 



1 1 

yo 1/ 



in 5(^+^2+), 
in ^(i?2+,i?2), 

in ^(^2,^2"), 



in 5(^2',^!), 

where 5'(7i,72) is the domain in Ob* bounded by 71 and 72 (see Figure H?7I) . 
Case V 

In this case (in comparison with the case IV) the branch point b* moves on the 0-sheet 
of 9^ (see Figure 12.111 and (12.651) ). As a result the arc E2 degenerates and after the 
normalization (14. 3 p of the initial Riemann-Hilbert problem (12.131) we have a new jump on 
51 instead of a jump on E2 (see Remark 12.321 -4). 



where 



on 51 C r^i ,0,2, 



$0 







1 






V 



which tends to the identity matrix as n ^ 00. In this case we open a global lens around 
Ai^2 in the domain bounded by and then, in the usual way, we open local lenses 
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around A2 and Ai. The model Riemann-Hilbert problem (the parametrix away from the 
branch points) for the function X for this case is identical with the corresponding problem 
of case I. Since b* is now on the 0-sheet of 9^, the local Riemann-Hilbert problem around 
b* is modified. Instead of the jump ( 14.50p in (]4.49p we have the jump 

J=}Wi, onAjnOfe*, , Sfe* := Ofo. n (A;^^^ U A* U5i). (4.52) 
(Dl, on 61 n Ob*, 

Thus the Airy function solution given in 04.511) for the problem (I4.49p . fl4.52p should be 
modified slightly: the non-trivial 2x2 block now is in the left-upper corner (instead of the 
right-lower corner in the case IV) and instead of ($1, $2, W1/W2) we now have ($05 
compare (I4.50p and (I4.52p . For the solution of an identical local 3x3 Riemann-Hilbert 
problem, see [SDl SHI iSl 13 EB [52] . 

Case III 

The proof for this case is just a repetition of the proof for the case V, with one simplifi- 
cation. Since Ai 2 = and Ai = AJ U 61 for this case, we do not need to open the global 
lens. Thus we are in the same situation as in the case I. The only difference is the solution 
of the local Riemann-Hilbert problem around 6*, which is the same as for the case V. 



5 Conclusion 

In this conclusion we highlight the main results of this paper. They are 

1. The classification of the sets A := {ai, fei; 02? ^2} such that the limiting counting 
measures for the poles and interpolation points (11.41) of the Hermite-Pade approxi- 
mants ([L2l) 

^{j) ._ D{i) p f . _ r)(j) 

for functions (12.41) 

fj G A{aj,aj;bj,Pj;n), 

are described by an algebraic function h of order 3 and genus 0, see (I1.19p . which 
give rise to measures A, /ii, and /i2 such that 

^Pn V2, z^^o) /^i, j = 1,2. (5.1) 



2. A universal vector-potential equilibrium problem (I2.83P for these limiting mea- 
sures. 

3. Strong asymptotics for the corresponding Hermite-Pade approximants. 
More precisely 
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In the definition of the classes of sets A we used a distinction in the formation of 
the system of curves 

r := {2 G C : Re jhj{z)dz = Re j hk{z) dz, j k, j, k = 0,1,2} (5.2) 

with a certain normahzation of the primitives. Since the genus of is the set F 
admits an algebraic parametrization. It is formed by the trajectories 



U 

1=1 



7£, ■= Zi{r]), i=l,...,6, 7] e[-2,2]. 



given by the branches of an algebraic function z{t]) of order six, when rj runs from 2 
to —2. The peculiarities of the behavior of these trajectories define different classes 
of A. Thus, starting from the input data A (i.e., the branch points ai, bi, a2, ^2) we 
have a finite number (because the genus of /i is 0) of algebraic functions h satisfying 

h^-?,^h + 2^ = 0, (5.3) 

with Il^{z) = {z — ai){z — bi){z — a2){z — 62)- From the coefficients of the equation 
(15. 3p we get explicit expressions of the coefficients of the equation for the function 
z{ri). Then, observing the behavior of the branches zlt]) when rj runs from 2 to 
—2, we can conclude to which class A belongs. Depending on the class we define 
(globally in C) the branches of the algebraic function h := {hQ,hi,h2} and the 
algebraic function 

$ := exp (^j dz^ = {"^o, "^i, '^'2}, 

i.e., we define the sheet structure of the Riemann surface 91 of the function h. 
Finally, the jumps of the function h on certain parts (depending on the geometrical 
case) of the contour F given in (15. 2p . give the densities of the limiting measures A 
and /Xi,/i2 in (15. ip . In particular 

dm = 7^ (/^o+(0 - ho-iO) d^, e e Ao C F, 
2m 

where Aq is a union of cuts which form the boundary of the domain of analyticity 
of the branch ho, i.e., Aq are the cuts of the 0-sheet fHo of the Riemann surface 9^. 

The vector-potential equilibrium problem (Theorem 12.3 ip does not depend on 
the geometrical class of A (universality). The equilibrium relations of this problem 
are considered on cuts Ai and A2 joining the pairs {aj,bj), j = 1,2. In addition, 
if Ai n A2 = Ai 2 7^ 0, an extra equilibrium relation is imposed on a curve E 
containing Ai 2, see (I2.83p . This equilibrium problem reduces to the known vector- 
potential problem for an Angelesco system when Ai 2 = 0- The cut Aj makes the 
function fj holomorphic. Thus the sets Ai, A2, Ai 2,-E are natural input data for 
the vector-potential problem and we expect that the limiting measures (15. ip of the 
Hermite-Pade approximants hold for a wider class of functions. We also recall that 
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we impose in this paper an extra analyticity condition on {fi, f^) when A12 7^ 0: 
we require that the ratio of the jumps of /i and /2 on Ai 2 

«(0 := y^— 7^(0, e e Ai,2 (5.4) 

J 2+ — J2- 

has a holomorphic (meromorphic) continuation from Ai 2 to the domain G 

u G H{G), dG = E. (5.5) 

The analyticity condition (15.41) gives a hnk of our vector-potential problem (12.831) 
with the known equilibrium problem for a Nikishin system. 

• For the derivation of the strong asymptotics for the multiple orthogonal polyno- 
mials Pn and for their functions of the second kind i?n (j = 1, 2), we use a 3 x 3 
matrix-valued Riemann-Hilbert problem. The increase of the order of the matrix 
functions (in comparison with 2x2 matrix-valued Riemann-Hilbert problems for 
the usual orthogonal polynomials) brings new features into the standard Riemann- 
Hilbert technology. One of these features is a new decomposition of the matrix jump 
on Ai 2 7^ which implies the opening of a global lens containing the domain G, see 
(15. 5p . This procedure introduces a new effect (in comparison with usual orthogonal 

polynomials) of oscillatory asymptotics on some curves C -E in C \ (Ai U A2) for 

(i) 

the functions of the second kind Rn and, as a result, there is an accumulation of 
zeros of Rn^ on E. 
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